Microscopic Theory of Vortex Dynamics in Homogeneous Superconductors 
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Vortex dynamics in a type-II superconductor is systematically investigated by the influence func- 
tional method. The irrelevant fermionic degrees of freedom are integrated out and their effects on 
the dynamics are treated in terms of the vortex coordinate. When an isolated vortex is moving 
against its background, forces proportional to the first order of vortex velocity on the vortex are 
calculated within the present formulation. The total transverse force on the moving vortex is ex- 
plicitly shown to be proportional to the superfluid number density and insensitive to impurities. Its 
equivalent expressions in terms of the Berry phase and the various summations of transitions be- 
tween quasiparticle (hole) states are discussed. At finite temperatures, due to the finite population 
of quasiparticle (hole) excitations above (below) the energy gap, there is a friction against vortex 
motion which diverges logarithmically in the low- frequency limit. Nonmagnetic impurities give rise 
to an additional friction from the core states which saturates to a value independent of the normal 
state resistivity in the dirty limit. In this limit, the coupling to the electromagnetic field does not 
change the conclusions if charge neutrality in the superconductor is maintained. Macroscopic con- 
straints on vortex dynamics by the second law of thermodynamics and by the fluctuation-dissipation 
theorems are also discussed. 



I. INTRODUCTION 

In a type-II superconductor, vortex motion is responsi- 
ble for a variety of low frequency transport phenomena. 
It is the only topological singularity whose dynamical 
properties are widely accessible to experimental studies 
in both classical and quantum -regimes, and its impor- 
tance has long been realizedE □ Despite decades of re- 
search, the theoretical agreement reached so far is very 
limited: At zero temperature, in the absence of any im- 
purity potentials, a vortex follows the local superfluid 
velocity. In the absence of a local superfluid velocity, 
when a vortex follows the motion of an external trapping 
potential, there is a momentum change in the superfluid 
transverse to the direction of vortex motion. In order 
to provide this momentum change, a force must be ap- 
plied by the external trapping potential to the superfluid 
through the vortex. The vortex experiences a transverse 
force proportional to the superfluid number density, bal- 
anced by the external force from the trapping potential. 
Beyond this simplest and idealized situation many as- 
pects of vortex dynamics have remained unsettled and 
even controversial. In the present paper we attempt to 
provide an influence functional formulation of vortex dy- 
namics from the microscopic Bardeen-Cooper-Schrieffer 
(BCS) theory, and a few detailed microscopic calculations 
under realistic conditions. 

The current microscopic understandings of vortex dy- 
namics in the presence of impurities and at finite tem- 
peratures may be classified into two different physical 
pictures, which are based on different theoretical ap- 



proaches and give contradictory results. In one picture 
the magnitude of the total force experienced by the vor- 
tex in the transverse direction is proportional to the su- 
perfluid number densityaQ The superfluid momentum 
change caused by the vortex motion is provided by an 
externally controlled trapping potential in the absence of 
a local superfluid. velocity, regardless of the existence of 
the normal fluid. eI It is an exact consequence of the global 
topological constraint on the vortex. The normal fluid at 
finite temperatures gives rise to friction for the vortex 
motion in the longitudinal direction. Furthermore, the 
global methods used in Refs. || and ^ indicate that the 
total transverse force is insensitive to random impurities, 
though there are additional frictional effects. In this pic- 
ture, in the absence of the externally controlled trapping 
potential, the pinning and friction should be used to ob- 
tain vortex motion perpendicular to the direction of an 
externally applied current. For the other picture, the 
essence of the results is that there are additional forces 
in the transverse direction of vortex velocity, provided by 
unbounded quasiparticle excitations or the normal fluid, 
by bounded vortex core stateSpj hji the substrate, or by 
a certain combination of themllllil The total transverse 
force is reduced, which is most clearly represented by 
the alleged gradual turning on the cancelation between 
two topological effects by a relaxation timefl the spectral 
flow of vortex core state transitions and the Berry phase 
counting far away from the core. To discuss this con- 
troversy from a detailed and straightforward approach is 
one of the main purposes of the present paper. 

We now state precisely the physical quantities which 



we are going to address. In the classical limit, we are 
looking for an effective equation of motion for a vortex. 
In two dimensions (2D), or for a straight vortex line in 
3D, the equation for a vortex specified by r^, takes the 
form of a Langcvin equation: 



„fi, = F(r„,t) - Br^ 



1] r,„ 



(1) 



Here rriy is the vortex mass, ij is the friction coefficient, 
B is the coefficient for the transverse force with z the 
direction perpendicular to the plane of vortex motion, 
and the fluctuating force f related to the friction force by 
the fluctuation-dissipation theorem. The force F contains 
all other forces which are not functions of vortex velocity: 
the force from the trapping potential, pinpings, the force 
due to an externally applied supercurrent,t3 the force due 
to other vortices, etc. We may classify the terms in the 
above equation into three types, according to the order 
of time derivatives of r„ : 

i). Forces contain no explicit dependence on any 
time derivative of vortex coordinate r.^, represented by 
F. These types of forces may be regarded as con- 
ceptually well understood, corresponding to the Born- 
Oppenheimer potentials, and are not controversial, 
though practically they can be difficult to evaluate. They 
are contained in our formulation, but will not be dis- 
cussed in the present paper. 

ii). Forces have a first-order time derivative of the vor- 
tex coordinate, the vortex velocity, represented by the 
transverse and longitudinal coefficients B and rj. Calcu- 
lating those forces explicitly is the focus of the present 
paper. It is our purpose to clarify the physical origins 
behind those forces, starting from a well-defined micro- 
scopic theory, the BCS theory, using a well-defined and. 
rigorous procedure, the influence functional method.cJ 
We will show in detail that the total transverse force is 
insensitive to details, and is proportional to the super- 
fluid number density, and present calculations leading to 
finite vortex friction contributions. 

iii). The term contains the second-order time deriva- 
tive of the vortex coordinate, the vortex mass rriy. This 
is also an uncleaf quantity, and the subject of the recent 
active study.Ljilj Though we believe our present formula- 
tion also provides the framework to address the dynam- 
ical eSects on the vortex mass, we will not explore them 
here.ll3 This term will be ignored by assigning the vortex 
acceleration r^, = 0. 

We organize the rest of the paper as follows. In Sec. II 
the total transverse force is studied from a macroscopic 
point of view. We first demonstrate from a thermo- 
dynamic consideration that the magnitude of the total 
transverse force should be proportional to the superfluid 
number density. A reduction from this value will lead to 
a violation of the second law of thermodynamics. Then 
we put the transverse force and friction in the context of 
fluctuation-dissipation theorems, and illustrate that the 
relaxation time approximation in the microscopic deriva- 
tions of vortex dynamics should be avoided. In Sec. Ill 



we first present a general formulation based on the BCS 
theory. Then we relate this formulation to that of the 
influence functional approach which has been proved to 
be rigorous and effective to calculate friction in quantum 
dissipative dynamics of a subsystem, where the total sys- 
tem is described by a Hamiltonian. A few general prop- 
erties of our formulation will be discussed. In Sec. IV we 
give detailed evaluations of both longitudinal and trans- 
verse correlations in the clean limit for arbitrary tem- 
peratures, which lead to both the friction and the total 
transverse force. In particular, detailed evaluations of the 
total transverse force from either extended state counting 
or core state transitions are given there, and are explicitly 
shown to be equivalent. In Sec. V the effects of impurity 
potentials are considered. We will show that the total 
transverse force will not be affected. However, impurity 
potentials strongly affect the core state spectrum, which 
leads to a contribution to the vortex friction in addition 
to that of extended states. In Sec. VI we show that the 
electromagnetic field does not affect the total transverse 
force and friction, under the condition that the charge 
neutrality in the superconductor is maintained. In Sec. 
VII some experimental tests are briefly discussed, and we 
summarize in Sec. VIII. 



II. THERMODYNAMICS AND STATISTICAL 
MECHANICS 

A. Force Balance and Thermodynamics 

The microscopic calculations which we will present 
later are unavoidably lengthy and technical. It may be 
helpful to obtain an overall picture and useful informa- 
tion (as much as we can) under general but elementary 
considerations. In this subsection we give a thermody- 
namic consideration of vortex dynamics to show that 
there is a constraint on the total transverse force, and 
in the next subsection the derivation of vortex dynam- 
ics will be put in the context of fluctuation-dissipation 
theorems. 

We may write down a possible equation of motion for a 
quantized vortex in the absence of impurities, taking into 
account the possible role of the normal fluid in the motion 
of the vortex, in the limit of the vortex acceleration equal 
to zero; 

hpsZ X (u„ - Vs) - D{uy - v„) - D'z X (u„ - v„) 



F*^^* = . 



(2) 



Here ps is the number density of the superfluid, u^,, v^ 
and v„ are the velocities of vortex, the superfluid, and 
the normal fluid with respect to the substrate or the wall 
of the container. Those velocities are independent vari- 
ables. The velocity dependences are only in first order in 
Eq. (2). The first term in the left-hand side of Eq. (2) is 



the Magnus forceja whose magnitude is proportional to 
the superfluid number density. The last term represents 
a possible external force on the vortex. The other two 
terms are possible contributions coining from the inter- 
action of the vortex with the normal fluid. Initially, both 
the normal fluid and superfluid velocities are set to zero. 

We will demonstrate that the conditions of force bal- 
ance and thermodynamics put a constraint on the value 
of D'. For this purpose let us imagine a torus-shaped 
tank filled with a superfluid, or a torus-shaped supercon- 
ductor film. The tank can be considered as a thermal 
reservoir to the superfluid. This implies that at finite 
temperatures there is also a normal fluid. After creating 
a vortex-antivortex pair, we keep the antivortex at rest 
and move the vortex to wind once with a small velocity 
Uy around one of the two circumferences of the torus, say 
Ly, in time ttotai before the annihilation with the antivor- 
tex. We take Uotai much longer than the relaxation time 
of the normal fluid such that the normal fluid velocity 
always stays close to zero, by transferring a possible mo- 
mentum gained from the vortex motion to the substrate, 
via the relaxation process represented by the normal fluid 
viscosity. A physical realization may be the electron- 
phonon interaction. Hence, the normal fluid velocity is 
always negligible comparing with the vortex velocity u^ 
which is an order of Uy « Ly /ttotai- As a result of the 
vortex motion, the momentum of circulating superfluid 
particles along the torus has been changed from zero to 
Psj = h/ Lx, because of the change of the winding num- 
ber of the superfluid. This occurs regardless how slow the 
vortex motion is. The kinetic energy of the superfluid 
has been changed from zero to E = psL^Lypl J2m* 
for a neutral superfluid, or when the effective magnetic 
screening length is larger than Ly for the superconductor 
fllm. Here m* is the effective mass of superfluid particles, 
and Lx and Ly are the circumferences of the torus. 

The total momentum change of the superfluid requires 
a force in the transverse direction of the vortex motion 

dP psLxLyPsj 
F^ = -7r= . — - = hpsUy , 

"I Hotal 

here P is the total momentum of the superfluid. Since 
the normal fluid velocity stays zero, no kinetic energy 
can be transferred from the normal fluid to the super- 
fluid. However, if the normal fluid would contribute to 
this force to superfluid by changing its internal energy, 
an additional transverse force on the vortex, —D'z x u„, 
arises. The magnitude of the external force in the trans- 
verse direction of vortex motion should be equal to the 
total transverse force according to Eq. (2), 

FT' = {hps - D')Uy . 

Now we are ready to consider the thermodynamic re- 
lations. The process of creating a vortex-antivortex pair 
and its annihilation after the vortex crossing one circum- 
ference Ly of the torus leaves only a flnite increase of 
superfluid circulation in the tank, corresponding to the 



change of winding number. The initial and flnal normal 
fluid velocities are zero. The increase of kinetic energy 
of the superfluid needs to be provided from somewhere. 
There are only two possible sources; the external trap- 
ping potential and the normal fluid. Here we need to 
be reminded of a signiflcant difference between the su- 
perfluid and the normal fluid: The superfluid carries no 
entropy, while the normal fluid does. Therefore, accord- 
ing to the second law of thermodynamics, EZl the super- 
fluid cannot gain kinetic energy by lowering the internal 
energy or entropy of the normal fluid. 

We need to consider the work performed on the system 
by the external force to move this vortex. In the longitu- 
dinal direction of vortex motion, the interaction between 
the normal fluid and the vortex gives rise to a vortex fric- 
tion —Duy. Thus the external force on the normal fluid 
in the longitudinal direction is Dutj. This friction does 
not dissipate energy. Rephrased alternatively, the energy 
dissipated can be arbitrarily small by taking the time to 
complete the process arbitrarily long, ttotai -^ oo. The 
process is then quasistatic. The normal fluid velocity is 
always negligible in the process because of its flnite vis- 
cosity. Thus the external force on the normal fluid in the 
longitudinal direction of vortex motion does not provide 
any work to the system. The external force acting on the 
superfluid will be able to provide enough work for the 
kinetic energy increase only if F^' > F±, i.e., D' < 0, 
which leads to the conclusion that the magnitude of total 
transverse force cannot be reduced from that determined 
by the superfluid number density. The work done by the 
external force is exactly equal to the kinetic energy in- 
crease if the magnitude of the external force in the trans- 
verse direction is the product of the superfluid number 
density, the Planck constant h, and the vortex velocity. 

The next question is whether or not the total trans- 
verse force on a moving vortex can be larger than that 
determined by the superfluid number density. If the nor- 
mal fluid would carry a vortex with a vorticity in the 
same direction as that of the superfluid, the answer to 
this question is positive. However, since we have assumed 
that the normal fluid is viscous, the vortex of the normal 
fluid will eventually disappear. This is true for a slow 
process whose time scale is much larger than the relax- 
ation time of the normal fluid assumed here. This con- 
sideration leads to that the total transverse force cannot 
be larger than the value determined by the superfluid 
number density. Combining with the thermodynamic ar- 
gument we conclude D' = 0. 

The above discussion has explicitly made use of the 
assumption of a flnite normal fluid viscosity. In case that 
the normal fluid relaxation time would be inflnite, that is, 
the normal fluid viscosity would be zero, a process which 
generates a vortex circulation in the superfluid would also 
generate a vortex circulation in the normal fluid. This 
would be the limiting situation of a dynamical process 
in which the internal relaxation time of the normal fluid 
is much shorter than its relaxation time to the substrate 
and the time scale for the process is between them. An 



example would be the creation of vortices by a magnetic 
flux in an ultraclean superconductor. In such a case, 
the normal fluid velocity will not relax to zero. Under 
this ideal condition —D' = hpn, corresponding to that 
the normal fluid has a vortex, which is what has been 
discussed in Ref. |l^. 

If impurities are present, a phenomenological equation 
of motion for the vortex may be written down if the im- 
purities are homogeneously distributed and vary only at 
a scale much smaller compare with the size of the vortex 
core. We have two more possible parameters from the 
vortex-impurity interaction: 



hpsZ X (u„ - Vs) - D{uy - v„) - D'z X (u„ - v„) 



duy — d'z X Uy 



= 



(3) 



Parallel to what we have discussed for the normal fluid 
case, d' must be zero in order for the external force to 
provide the energy gain needed by the superfluid. The 
impurities cannot provide energy to the superfluid either 
by lowering their internal energy or entropy because of 
the second law of thermodynamics. We note that impu- 
rities introduce another contribution to the normal fluid 
viscosity. |-||-| 

The microscopic global considerationsoQ have already 
suggested D' = and d' = 0. The conclusion here will 
be borne out by detailed and independent microscopic 
calculations in the following sections. 



B. Friction and Fluctuation-Dissipation Theorems 

The derivation of the equation of motion for the vortex 
is different from the usual linear-response theory. In the 
linear-response theory, a driving force is given, that is, 
the Hamiltonian is known, and we look for the average 
responding velocity. It is a calculation of conductivity or 
mobility. In the present case the force on a moving vortex 
is the unknown quantity which we need to find out. The 
vortex velocity is, however, readily defined through the 
vortex coordinate. It is a calculation of resistivity or fric- 
tion. To appreciate this difference, we will examine the 
different correlation functions involved in these different 
types of calculations and their relationships. The focus 
point in this subsection in on the condition for using the 
relaxation time approximation. 

When the normal fiuid is at rest, the vortex motion is 
governed by a classical Langevin equation with param- 
eters to be determined microscopically. This equation 
has the same form as a classical electron moving in a 
magnetic field. We adopt the language in transport the- 
ory i|a,make it easier to relate to the early work in that 
fieldMM 

We start by considering a classical charged particle in a 
magnetic field obeying a generalized Langevin equation: 



mu,{t) = dt' ry,fe(i - t')uk{t') + Fr\t) 

Jto 

-Be,kUkit) + Mt) . (4) 

Here i = x or y, u(t) = {ux{t),Uy{t)) is the velocity of 
the particle, m is the mass, F^^*(i) = {F^''\t), F^''\t)) is 
an external force, f(t) = {fx{t), fy{t)) is a random force 
which simulates the effect of the thermal reservoir. The 
Einstein convention of the repeated indices as summa- 
tion has been used. BeikUkit) represents the transverse 
force, the Lorentz force —u{t) x B in the Langevin equa- 
tion with the magnetic field taken along the z-direction. 
The matrix ri{t — t') = {rjij} represents friction in both 
longitudinal and transverse directions of the particle mo- 
tion. Its possible finite off-diagonal elements will change 
the effect of the original Lorentz force on the particle. In 
addition, we have 



< m) > = , 

<u,(to)/,(to + t) > =0, t>0 



(5) 



< u^{to)uj{to) > 



knT. 



The first equation is obvious: no average fluctuating 
force. The second one is due to the causality and the 
last one is due to the equipartition theorem. 

If the Hamiltonian of the particle is known, the prob- 
lem of particle responding to a perturbation can be for- 
mulated in two different but equivalent ways. We can 
calculate the velocity of the particle while the applied 
force is given. In such a case, it is to obtain a conduc- 
tivity or mobility formula. The conductivity or mobility 
may be obtained by the Nakano-Kubo's formula, a cal- 
culation of velocity- velocity correlation function. It may 
also be obtained by solving thaiBoltzmann equation in 
the presence of an electric field.u Otherwise, we can con- 
sider a given velocity for the particle and calculate the 
applied force needed to maintain this motion. It is to 
obtain a resistivity or friction formula, i.e., calculating 
an electric field needed to maintain the given current. 
The derivation of vortex dynamics belongs to the second 
kind, where we consider a steady motion of the vortex 
and calculate the external force acted on the vortex. Un- 
fortunately, we do not have the choice to formulate vortex 
dynamics in superconductors in terms of conductivity or 
mobility formula because the effective vortex Hamilto- 
nian is unknown a priori. 

Introducing a Laplace transform 



VM 



dt 



*#) 



the mobility is given from Eq. (4) in the limit to 

by 

IJ.[uj] = {imuj + rj[uj] + iayB)^^ . 
Here the mobility /i[a;] is defined through 



(6) 



< Ui[uj] >= /iy[tj]r^^ [uj] 



with an applied external force F (i) = F [tjje'"*. 

Defining the velocity-velocity correlation function ma- 
trix 

with Vlij{t = 0) = 5ijkBT/m according to Eq. (5), 
the mobility is related to the velocity- velocity correlation 
function 



fl[UJ\ 






(7) 



This is the 'first! fluctuation-dissipation theorem de- 
scribed by Kubo,E£l equivalent to the Nakano-Kubo's for- 
mula for the electrical conductivity. 

It is easy to demonstrate that the relaxation time ap- 
proximation can be valid in the Nakano-Kubo's formula. 
Without the thermal reservoir, the velocity-velocity cor- 
relation is given by 



U[uj] = {imoj + layB) ^m U{0) . 



(8) 



When using a relaxation time approximation by the stan- 
dard rule, iu) —f iu> -\- r][uj]/m and substituting it into Eq. 
(8), we find the velocity- velocity correlation under the 
relaxation time approximation is given by 

U[ijj] = {imuj + ri[Lj] + iayB)^^kBT , 

which is exactly the same as the one obtained by the 
rigorous calculation, Eqs. (6) and (7). Therefore, the 
relaxation time approximation can be a valid one for 
velocity- velocity correlations when used in a conductivity 
or mobility formula. 

The resistivity or friction formula is known tjX|hfj_dif=. 
ficult and it is worthwhile to examine it closely .oEilo'LJ 
First, we calculate the total force-force correlation func- 
tion matrix 

Tij{t) = m^ < Mj(io + t)uj{to) > . 

Taking the Laplace transform, using the translational in- 
variance in time 

< Uj(to +t)uj{to) >= - < ^^(to + t)uj{tn) > , 

and the total force-velocity correlation function 

m < Ui(to + t)uj{to) > [lu] = -~mLlij{Q) + iraujUij[ijj\ , 

we have 



TH 



-iBa,, 



imoj - 



'7N 



iayB 



knT. 



(9) 



rn 



-iBay ksT , 



(10) 



and is independent of the frictional coefficient 77. 

The random force-force correlation matrix is defined 



as 



T^v{t)^<hito + t)fj{to) > 



(11) 



^From the Langevin equation, Eq. (4), we can express 
TZ{t) in terms of total force-force, total force-velocity, 
and velocity- velocity correlation functions. Taking the 
Laplace transform and integrating by part, we obtain 



n[Lu] = f][u] mU{0) = ri[uj] kgT , 



(12) 



or rj{t) — TZ{t)/{kBT). This is the 'secpad' fluctuation- 
dissipation theorem described by Kubollj We emphasize 
that the generalized frictional coefficient ri{t) is directly 
given by the random force-force correlation. The fric- 
tional coefficient matrix r](t) has no off-diagonal part if 
the random force- force correlation matrix has not. This 
fluctuation-dissipation theorem allows us to obtain some 
general properties of the generalized friction. For exam- 
ple, for a charged particle described by a single relaxation 
time in the Boltzmann equation moving in a magnetic 
fleld, there will be no frictional effect on the force in the 
transverse direction. 

Next we consider that the particle is moving at a given 
velocity u{t) and flnd out what is the external force 
needed to sustain such a motion. This is exactly the 
situation which we encounter in vortex dynamics and it 
is equivalent to the calculation of resistivity or friction. 
^From Eq. (4), the average force is given by 



<Fr'\Lu] >= 



ri[uj] + i(7yB)i 



(13) 



which is trivially identical to the reciprocal of conductiv- 
ity formula. Obviously, this process does not provide us 
an independent way of calculating resistivity. 

However, if we are only interested in the average force 
< F*^^* > in the DC limit, we do have an alternative 
resistivity formula. After taking w ^ and using Eqs. 
(10) and (12), Eq. (13) gives 

< F^ > [0] = fc^ (7^^.- [0] - ^^0 [0] ) u, [0] . (14) 

Taking 77 to be a scalar, the external force can be written 
into a more suggestive form. 



< F'^^'p] >= 77[0] u[0] + u[0] X B , 



(15) 



In the limit 



0, it is reduced to 



where the longitudinal component depends on TZ[Q\, the 
random force-force correlation function, and the trans- 
verse component only on T[Q\, the total force- force cor- 
relation function. Equation (14) is the DC resistivity 
formula. It provides a direct way to obtain DC resis- 
tivity from force correlation functions. The straightfor- 
ward interpretation of Eq. (15) is the force-balance: The 
externally applied force to keep the constant velocity is 



equal in magnitude but opposite in sign to the sum of 
the frictional and the Lorentz forces. 

We will show that the relaxation time approximation is 
invalid when used in the force calculation. We start with 
the force-force correlations. Without thermal reservoir, 
the random force correlation is zero, that is, TZ{t) = 0. 
If we switch on the effect of a thermal reservoir by using 
a relaxation time approximation icu —^ ito + ri[u!]/m, the 
random force correlation is still incorrectly set to zero. 
This shows that the relaxation time approximation can- 
not be used to calculate the random force correlation. 

The total force correlation without thermal reservoir 



IS 



Tij [Uj] = -iBCTy 



(muj)^ 



imuj + iayB y ^^ 



mUkjiO). 



(16) 

When we switch on the thermal reservoir by using a re- 
laxation time approximation ilu ^ ilu + 7y[w]/?n in Eq. 
(16), we have 

.FM = -iB<jy + imu + 77M - . \ / !^ keT 

\ tmu) + r][uj\ + layB J 

This is a rather complicated expression. In the limit 
Lu << riluj], or ujT << 1, we can simplify it to 



T[0] = 



B 



1 + (c^0t)2 



{ujqt - iay{ujQT)^) ksT , (17) 



with wq = B jm. Here t = ?7i/77[0] is a relaxation time. 
Let us use the resistivity formula Eq. (14) to calculate 
the external force needed to keep the particle moving 
with a given velocity. With 'R,\u)\ =0 and JF[0] given by 
Eq. (17), the external force is 



< F'=^*[0] >= - 



CJqT 



(wot) 



r(-Bu[0]-hworu[0] xB). 



(18) 



These results have no connection at all to the rigorous 
results shown in Eq. (15). Even the sign for the lon- 
gitudinal force is wrong. Evidently, the relaxation time 
approximation cannot be valid in such a calculation, be- 
cause such an approximation in the force balance equa- 
tion cannot consider the random force properly, and leads 
to results violating the fluctuation-dissipation theorems. 

By the simple and exactly solvable model, we have 
demonstrated the essential conditions for the validity of 
the relaxation time approximation in velocity-velocity 
correlation function calculations, and for its invalidity 
in force-force correlation function calculations. We refer 
readers to Refs. [l9| , pi| , p(]| ,p2|for more sophisticated discus- 
sions in the context of the Green's function or Boltzmann 
equation. 

With a redefinition of constants B ~ hps with h be- 
ing the Planck constant, ps being the superfluid parti- 
cle number density, and ujq = eo as the core-level spac- 
ing, this force becomes the same in magnitude as the 



one which appeared in the derivation of WrtPii dynamics 
using the relaxatioit-time approximation^ B'Ell including 
the same sign error .t3 The inappropriate use of the relax- 
ation time approximation in vortex dynamics in d-jwave 
superconductors has also been pointed out recentlyBj In 
the following we show how to obtain the vortex friction 
without the relaxation time approximation, and demon- 
strate at the same time that the total transverse force 
remains unchanged as dictated by the topology. 



III. VORTEX DYNAMICS IN HOMOGENEOUS 
BCS SUPERFLUID 

A. Formulation of the Problem 

We present now our microscopic derivation, from the 
standard BCS Lagrangian for s-wave pairing in the imag- 
inary time path-integral formulation of the influence 
functional method. The connection of the total trans- 
verse force to the Berry phase is straightforward in this 
formulation. We believe the present formulation has 
some advantages: a transparent crossover from the quan- 
tum to the classical description via the semiclassical ap- 
proximation, and a flexible treatment of the general dissi- 
pative effect arising from the integration out of irrelevant 
degrees of freedom, fermionic quasiparticles, and holes. 
The relevant degree of freedom is the vortex coordinate.tj 

We consider a neutral fermionic superfluid first. The 
coupling to electromagnetic fields will be discussed later. 
The Lagrangian is given by 



Lbcs = X! ^'^(^' ^'^ ( ^^^ ~ ^^ 



2to 



V{x) 



-5-0I {x, r)-0| (x, T)V'i {x, r)V'T {x, t) , (19) 

where ip^ describes electrons with spin a = (t,i), Pf is 
the chemical potential determined by the electron num- 
ber density, V{x) is the impurity potential, t/o is the 
trapping potential, and x = {x, y, z). A vortex at Xy has 
been assumed through the trapping potential. A more 
explicit implementation of the vortex coordinate will be 
discussed after Eq. (39). The partition function is 

Z = I?{a;^,,V'^'0}exp< -- dr cPxLbcs \ , 

(20) 

with (3 = l/ksT , and (Px — dxdydz. Inserting the 
identity in the functional space, 

1= /"l?{A*,A}expJ-| / dr I d^x 



V'iV'T + -Mx,t) 



into Eq. (20) we have 



Z= /"I?{x„,V^V',A^A}expi-i / dr f d^x 



{i;\,i;i){ndT + H) 



h/3 



dr I d^xlA^^ 




1 

^g Jo 

Here the Hamiltonian is defined as 
H(A,A*)=(f, _^, 



(21) 



(22) 



with H = -{n'^/2m)V'^ - fip + V{x) + Uo{x-x^). 

Exactly integrating out the electron fields V'J and 4'^ 
first, then integrating out the auxiliary (pair) fields A 
under the mean-field approximation, one obtains the par- 
tition function for the vortex 



Z = Vlxy} exp < - 



S, 



eff 



(23) 



with the effective action 



Seff 



1 



fi/3 



= -Tr In G-i + — / dr / d^xjAp , (24) 
n- rig Jo J 

where Tr includes internal and space-time indices, and 
the Nambu-Gor'kov Green's function G is defined by 

[Tidr + n)G{x, t; x\t') = S{t - t')6^ {x ~ x') , (25) 

together with the BCS gap equation, or the self- 
consistent equation. 



A{x,t) =hg Gi2 (x, r; x, r) 



(26) 



In the presence of impurity potentials, the averaging over 
them is implied in Eq. (26), unless explicitly specified. 

Since the effective action is a smooth function of vortex 
coordinate in the functional space of {xv{t)}, we consider 
that the vortex has made a small move from its original 
place xo, which allows a small parameter expansion in 
terms of the difference between the vortex position x^ 
and xq. We look for the long-time behavior of vortex 
dynamics under this small parameter expansion. In the 
final step, the forces on vortex are to be calculated by 
varying Lagrangian to this small motion. The influence of 
the eliminated degrees of freedom on the vortex dynamics 
will then be obtained. As an example, for the mean- 
field value of the order parameter, this small parameter 
expansion to the second order is 



A{x,T,Xy) = [ l + SXy{T) • Vj,o +-i6Xv{T) ■Vxo)'^ ] X 



Ao(x,a;o) 



(27) 



Here Sx^ = Xy -- xq- In Eq. (27) we have used the fact 
that when Xy — xq the vortex is static. The effective 
action for the vortex to the same order is, after dropping 
a constant term, 



Seff 1 . 



'e// 



i-hp 



SxyVx^A*SxyVx^Ao (28) 



with 



E' = Sxy ■ V 



V ^ XQ \ A * 



Uo Ao 
A5 -Uo 



Sxy ■ Vx.Ho ■ (29) 



Here the Hamiltonian Tig = Wl^j^^ajf, for the static vortex 
at Xq, Gq is the Nambu-Gor'kov Green's function with H 
replaced by Tig, the gradient V^o is with respect to xq, 
and G-i = Go ^ + S'. 

Now we construct the Nambu-Gor'koy Green's func- 
tion Go following the usual procedure. c3 First, we con- 
sider the eigenfunctions of Ho- The stationary equation, 
the Bogoliubov-de Gennes equation, is 



Wo*a(2:) ^Ea-^aix) , 



(30) 



with 



*a(x) 



Ua{x) 
Va{x) 



No confusion with the vortex velocity in Sec. II should 
arise here. 

Given the eigenfunctions of Eq. (30), Go can be ex- 
pressed as 

G„(x,r;.'.,').g^g;;:^*.W*.(..'M31) 

Here LUn = ni: /T%j3, with n odd integers. 

Direct substituting of Eq. (31) into Eq. (28) leads to 



S, 



eff 



1 /• -iw„(r-T') -i(^„,(T'-T) 

' d'xd^'x'dTdT' Y. ^-- ^X 






' 2{hpf 

■^i{x')5xy{T') ■ Vx,no{x')-^o.'{x') X 

*l,(a;)fe„(T) • VxQno{x)^a{x) 

1 /■'''^ /■ 
+ ^ / dT d^x6xy{T) ■ V^oA* Sxy{T) ■ V:,„Ao . (32) 



Keeping only terms relevant to vortex dynamics and as- 
suming global rotational symmetry after summing over 
all the states, we have 



'e// 



1 /■ , Q \-^ g— Jw„(r-T') g-'iw„,(T'-r) 

(PxcPx'dTdr' > ^ = : :^ X 



■^ Z?ICJ„ — Ea iflijJn' — Ea' 
na,n'a' 



2{h(3Y 

5x,{t') ■ 6x,{t) + (*l(x')V,„Ho*a'(a;')x 
*l'(a;)V:roHo*a(x)) • z (fa„(r') X fa„(r)) • z 

With a rearrangement, finally we arrive at 



R ~ c», one may consider the scattering states. In 
this case the Dirac 6 function normalization for extended 
states should be the better choice. Furthermore, Eq. (30) 
has the property that if 



(33) 



then 



Ho*(a;) = E ^{x) , *(a;) 



Ho^ix) = -E *(a;) 



v*ix) 
-u*{x) 



(37) 



5, 



eff 




l\\2 



dT'F||(r- T')\Sxy{T) - Sxv{t')\ 



dT'Fj_{T - T'){Sxy{T) X dx^ir')) -z} , (34) 



There is no specific assumption about the Haniiltonian 
H{H*) in Eq. (22) for this identity. There is another 
important property implied by Eq. (30). Since both the 
Hamiltonian Tio and its eigenfunctions are the function 
of the vortex coordinate at xo , taking the derivative with 
respect to xq at both sides of Eq. (30), we have 



with 

F|| (r-r') 
1 



2hl3'^ 



d^xd^x' Y^ 



g-icj„(r-r') g-iw„/(r'-r) 
iTlLOn — E„ itkUn' — En' 



*l(a;')V,„Ho*a'(a;') • *|,,(x)V,„Ho*a(a:) , (35) 

and 

^ ^ / d^xd^:,' y 1- 1 X 

W J ^^ ihiVn — Ea iriUJri' ~ -^a' 

na.,n'a.' 

(*l(a;')V,,Ho*a'(x') X *^,(a;)V,„Ho*„(x)) • z. (36) 

Equation (34) has the form of riafluence functional 
in quantum dissipative dynamics .Ii3 Please note that 
/„ i^l|(r — T')dT = 0. Therefore, there is no so-called 
'counter-term' in Eq. (33) as discussed in Ref. |l^. Here 
we have generalized the influence functional to include 
the transverse force as a response from the environment. 
Before proceeding to evaluate these correlations, we 
discuss some properties of the wave functions of the 
Bogoliubov-de Gennes equation, which will be used later. 
First, because Ti.Q is Hermitian, all its eigenstates form a 
complete and orthonormal set, that is. 



d^X'^i{x)-^a'ix) ^Sa,a' 



and 



^vl/„(x)*l(:r') 



Here ^^(a:) = {u* (x) , v* (x)) and the wave function ^(a;) 
is normalized to 1 over a cylinder of radius R and length 
L, the box normalization. In the thermodynamic limit, 



(V,oH)|*„. 



n\w,„'$^,)^E^,\w,,^^, 



Multiplying both sides of this equation by (^q|, and us- 
ing the relation that (^q,|7Y = E^i'i'al, the Hermitian 
conjugation of Eq. (30), for a 7^ a' we have 



d''x-^i{x){W.,,no)'^a'ix) 

- {E^, - E^) j d^x-^l{x)V,„-^^,{x) , (38) 



with 



(V,„7io) = V,, 



Uo A 
« > A* -C/o 



Here we have used \7^^Ea' = to get Eq. (38), un- 
der the assumption that the system is homogeneous. 
Hence there is no the vortex velocity independent po- 
tential for the vortex arising from Eq. (38), that is, 
no Born-Oppcnheimer type potential, in accordance with 
the present purpose of looking for the effects which are 
first order in vortex velocity. Starting from the Hermi- 
tian conjugate of Eq. (30), taking the derivative with 
respect to the vortex coordinate we have 

Then multiplying this equation by \'^a') we have 



d^X^Ux){W.,„Ho)'^a'{x) 

-{E^, ^ E^) j d^xV^,^i{x)^o.'{x) . (39) 



We note that both Eqs. (38) and (39) are exact, follow- 
ing from the general property of Eq. (30). They relate 
the transition elements of the Hamiltonian after the dif- 
ferentiation with respect to a parameter to the connec- 
tions between wave functions. Though the wave func- 
tions have to be determined as an eigenvalue problem. 



the usefulness of Eqs. (38) and (39) is that it allows one 
to concentrate on wave functions instead of the original 
Hamiltonian, which is particularly convenient in the dis- 
cussion of certain topological properties described better 
by wave functions, such as a vortex in a BCS superfluid 
here. In the rest of the paper, we will take the trapping 
potential to be zero, C/q ^ unless specified, and de- 
termine the vortex position self-consistently through the 
gap equation, Eq. (26). 

For the convenience of calculation, sometimes we wish 
to use V^a(a;) instead of Si'xo'^aix) in the expression. 
It can be done in the following way. We split the gap 
function, or the order parameter A into 

A = A(a;-xo) + A'(x,a;o) , 

where A is a smooth part of the self-consistent potential, 
A' is the fluctuating part for a given impurity configu- 
ration. The impurity average gives < A' >= 0. In the 
presence of impurity potentials the gap function A may 
differ from the one in the clean limit. The Hamiltonian 
becomes Ho =7^0 + STi, with 



Ho — 



Ho A 
A* -m 



where Hq 



-{li/2m)V^ - liF: and 

on - I ^,, _yf^^^ 



(40) 



(41) 



Using 



d''x^i{x)V{Ho^„,{x)) 



(fx [^i{x){VHo)'^c.'{x) + *i(x)7^oV*„-(a;)] , 
and (V + y xq)H.o = 0, and defining 

we have the desired relations 

<fx^i{x)Vx„^a..{x) 



and 



{Ea' — Ea 

= -iE„,-E^) I (fx^l{x)V^^,{x) 
+ I d'x-^i{x)VH'-^c.'{x) 

{E^,~E^) f d^xV^,^i{x)^c.'{x) 
= -{E^,-E„) I (fxV^l{x)-^^,{x) 
(f'x-^l{x)\/H'-^a.{x). 



(42) 



(43) 



The last part is obtained from Eq. (42) by a partial inte- 
gration, which can be carried through because the wave 
function is normalizable, either by the box normalization 
or by the Dirac delta function. 



B. Longitudinal Correlation 

We now discuss the general properties of the longitu- 
dinal correlation function, Eq. (35). We find 



E 



n,n' 



g-icj„(r-T') g-lCJ„/(T'-T) 
itiUJn — Ea ihuJn' — £■«' 

g-'i(w„-W„/)(T-T') 



Et 



', iiJiUJn - ?iW«') - {Ea - Ea') 



1 



1 



iflUJn' ~ Ea' iflUJn — Ea 



= Y.Pifo.'-uh 

^ ^ 7 



-'^,^'){T-T') 



i{huJn - huJn') - {Ea - Ea') 



(44) 



after using 



E' 



-. — iuJri^ 



lhcj„ ~Ea ] -(3{l-fa), 6 = 0+ 



with the Fermi distribution function /„ = 1/(1 + e^^"). 
To complete the calculation, we also need 



E 



C0s[(^n - UJn'){T - t')] 
\ i{huJn — fujJn') — {Ea ~ Ea' 



o cosh 


h 


^(M-|.- 


-.'!)] 


2 


sinh 


{Ea~Ea')§ 





+ :S{t-t): . (45) 



Here : 6{t) : is a periodic delta function with period 
hf3. The term with 

J2n-n' Sin[(Wri " W„')(r - T')]/[i{huJri " huJn') " {Ea - Ea')] 

is zero inside the double imaginary time integration in 
Eq. (34), because the integrand is an odd function of 
T — t' . Dropping the periodic 5 function, whose contri- 
bution is zero in Eq. (34), we are ready then to write 
down the longitudinal correlation function as 



cosh 



F||(t) == - / dujj{uj) 



UJ 



(f-M) 



sinh 



hf} 



(46) 



with the spectral function 

j{^) - 7 E'^(^^ - 1^" - ^"'1)1/"' - /" 

d^x^l{x)\/x„Ho'i'c.'{3 



(47) 



It is interesting to point out that in terms of the spectral 
function the longitudinal correlation function, Eq. (46), 
is in exactly the same form of the influence functional in 
quantum dissipative dynamics .E3 The apparent difference 
is that the spectral function in Rcf. n3 has been obtained 



by integrating out a set of independent harmonic oscilla- 
tors (bosons), while here it has come from the elimina- 
tion of independent fermionic modes determined by the 
Bogoliubov-de Gennes equation. 

It is important to remember that in order to have a 
smooth spectral function J (to), the thermodynamic limit 
must be taken first before the implementation of the 
S function in Eq. (47) £3 Otherwise the spectral func- 
tion consists of sum of discrete delta functions, and there 
would be no dissipation. This limit procedure is in ac- 
cordance with the requirement in nonequilibrium statisti- 
cal mechanics: The thermodynamic limit must be taken 
first to have well-defined low-lying modes in the zero- 
frequency limit. After this consideration of the thermo- 
dynamic limit, in the low- frequency limit the spectral 
function may have the following generic form: 



J{lu) =770;" 



0^ 



(48) 



with s > 1 being the super-ohmic case, 1 > s > being 
the sub-ohmic case, and s — 1 being the Ohmic case, fol- 
lowing from the influence functional formulation of quan- 
tum dissipative dynamics .E3 For the physically important 
Ohmic case, the longitudinal force, friction, is given by 
— ?7Vy, and from Eq. (47) we have the frictional coeffi- 
cient 






(49) 



This equation is the familiar Fermi Golden rule for dissi- 
pation. The matrix elements of ^ xqT^o are well behaved. 
If we use Eqs. (38) and (39) to re-express the frictional 
coefficient 77 in terms of the overlap integral between the 
wave functions \Wxo^a) and |v1/q'), we turn it into the 
form of ratio 0/0 when E^ — E^i -^ 0. Then attention 
should be paid to the divergence of the overlap integral 
when \Ea ~ Ea'\ -^ 0^. This limiting behavior has been 
discussed in Ref. ^ and we refer the reader to the Ap- 
pendix for a detailed discussion. Equation (49) clearly 
shows that the coefficient of friction 77 is determined by 
low-energy excitations such as phonons, extended quasi- 
particles, and bounded core quasiparticles when their en- 
ergy spectrum is smeared out by impurities. The equiv- 
alence of Eq. (49) in the context of vortex dynamics 
to a more conventional partial wave phase-shift analy- 
sis has been discussed in Ref. g^ for a few well-defined 
situations. A more formal discussion can be found in 
Ref. ^. It may also be instructive to mention here that 
the friction experienced by a moving object in a normal 
Fermi liquid has been analyzed in the infiuence functional 
approach. c2l Those considerations suggests that nonzero 
extended states friction contributions exist, as will be 
borne out in detail in the next section. Finally, it should 
also be pointed out that Eq. (49) is a special case of Eq. 
(47). It will not pick up any super-Ohmic contributions, 
and will give infinity for any sub-Ohmic contributions. 



If such cases occur, we need to return to the general ex- 
pressions, Eqs. (46) and (47). 

To close this subsection, there are two general remarks 
in order. First, the present result of expressing the fric- 
tional coefficient in terms of low-lying excitations is in 
accordance with Landau's quasiparticle picture: In the 
zero-frequency limit the lifetime of those excitations ap- 
proaches infinity. Those excitations give an exact de- 
scription of the dynamics of the whole system in this 
limit. Secondly, to relate to the discussion in Sec. II. B, 
the spectral function given by Eq. (47) completely deter- 
mines the spectral representation of the second kind of 
fiuctuation-dissipation theorems. Equation (47) is indeed 
a quantitative description of dissipation. 



C. Transverse Correlation 

To obtain physically more transparent expressions for 
the transverse correlation function in Eq. (36), we use 
Eqs. (38) and (39) to rewrite it as 






^(t-t') -i(^„/(r'-r) 



hB"^ J ^-^ ihjJn — Ea iflUJni — E^' 

na.n' a' 

{E^-E^,fx 
{^i{x')Vx,^o.'{x')) X (Vx,^l{x) ^c.{x))] -z. (50) 



According to Eq. (44) and (45), 



2^ in 



-ii^ujT^—ijj^i '){'t—t'') 



\ {ihuJn - Ea){ihuJn' - Ea') 
g-iil)„ (r-r') 



n 

2 

cosh 



ihoJn - {Ea - Ea') 



iU' - /a) 



-1 + 
Ec-E, 



-dr^T' 



'{!f-\r-r'\ 



sinh 



{E^-E^,)fi 
2 



(/a' - /a) . 



Because of the symmetry with the interchange of a and 
a', the —1 term in the above square brackets does not 
contribute to the transverse correlation function, and we 
have 

F^{t~t') 



= - / SxCpx' ^ dr-T' 

a,a' 

{Ea ~ Ea')ifa' - U) ^ 



cosh 


\E^-E, 


^(f-k 


-'!)] 




sinh 


"(B„-_E„,)/3l 
2 





:*i(a:')V.o*a'(x')) X (\7,^^>l{x) ^a{x) 



z . 
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The corresponding term in the effective action, Eq. (34), 
is then 



up i-np 

dr / dT'F±{T - T'){SXy{T) X SXy{T')) ■ z 

Jo 



^ dT dr' d^'xd^'x' Y, 



cosh 


\E^~E, 

h 


-(M-|r 


-'!)] 




sinh 


2 





dr-T' 

\^i{x')\/,,^^,{x')) X (V,„<,(x) *„(a;) 
((5x„(r) X dxv{T')) ■ z . 

= -^ dT dr' d^xd^x'J2 

Ea — E, 



dr-T' CXp I - 



"\-T'\\sgn{E^-E^, 



{^l{x')V,,^^,{x')) X (V^,^l,{x) ^^{x) 
((5xi,(r) X 5xv(t')) ■ z . 



z X 
(51) 



In the last equahty we have used the periodicity of the 
function 6xy (r) = 5xv {hp + r) to turn the hyperbohc 
function into the exponential function. Now we look for 
the slow motion expansion to the leading order in ve- 
locity: Sxv{t') — Sxv{t) + Sxv{t){t' — t). Substituting 
this expansion into Eq. (51), after the integration over 
(r' — r) we have 



1 



hp 



dr 



dT'F±{T - T'){SXy{T) X 6Xy{T')) ■ z 



hp 



dT B [5xv{t) X 5xv{t)] ■ z , 



with the quantity B which determines the transverse 
force defined as 

B = i\Y.^U'-U)jd'xjd^x' 

z-(*^(a:')V,„*„.(rE')x V,„*J,,(x)*,(x)) . (52) 

We demonstrate next that the contribution to the 
transverse correlation can be evaluated by counting ex- 
tended states contributions. First, we regroup terms in 
Eq. (52): 

5 = *^-^-EH2E«- 



2 Z-^ ^ ' Z-^ R ihuJn ~ Ea 

a,a' n 

d'^x j d''x'^a{x)^l{x')V^,^a,{x') X V^,^l,{x) 



,^2-trY,i-)^J2^ 



a n 

d^x{'sj^„^a{x) X V:,„<(X)) , 



(53) 



because that {^a'} form a completed set. Here tr stands 
for summing over spinor indices. The replacement of fa 
by the summation is to take care of the delicate equal- 
time limit in the trace: 5 — Q^ for spin up and i5 = 0+ 
for spin down in Nambu spin space. We encounter such 
a choice of time limit only in the case of taking the trace 
of the Nambu-Gor'kov Green's function directly. This 
choice will not be there if we only need to take trace of 
higher powers of the Nambu-Gor'kov Green's function, 
e.g., trGQ. This implies that the functions containing 
occupation numbers {fa} in their differences are well 
defined. Therefore, Eq. (52) can be safely used if we 
directly substitute in eigenstates of the Bogoliubov-de 
Gennes equation. An alternative natural way of deriv- 
ing Eq. (53) is to leave the summation over cOn in place 
throughout Eq. (50) to Eq. (53). 

After substituting Eq. (53) into Eq. (52), we will write 
it explicitly in terms of the eigenstates of Bogoliubov-de 
Gennes equation. Since 



1 „ — iuJnS 



d^xV^g^aix) X y^o'^iix) 



13 ihiOn ~ E, 

= E / d^x{fa'^XQUa(,x) X V XQU*a{x) 

-(1 - fa)^xoVaix) X V^o<(a;)} , 
we obtain 

B = -ihz ■ E / d^x{faVxoUa{x) X VxoUa{x) 
a '' 
-(1 - faWxoVaix) X ^xoVaix)) . (54) 

After using Wx„ -^ —V, we ©ra,luate Eq. (54) with the 
help of the current definitionj29 

J = ^ X! {fa'U*a'^Ua + (1 - /a)WaVi;*} + C.C. . 



Equation (54) becomes 



B 



dx z ■ {V X j) 



dl-j 



X — Xq I — >-00 

2TTnps{T) . 



(55) 



Here we have used the fact that the current is zero at 
the vortex position, as the explicit calculation in Ref. pG 
has shown. It can be understood as the requirement of 
quantum mechanics: At the phase singular point, the am- 
plitude of any wave function carrying this singular phase 
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must be zero. In reaching Eq. (55) we have also made 
the assumption that the vortex does not have a normal 
fluid circulation. The normal fluid is in equilibrium with 
the substrate or the walls of the container. 

/^From Eq. (55) one may conclude that when counting 
the contribution from individual states, only extended 
states give rise to the contribution to the transverse re- 
sponse, because the loop of the line integral can be cho- 
sen arbitrarily large to make the core state contributions 
arbitrarily small. It corresponds to the fact that only 
extcndejd states can contribute to the Berry phase of the 
vortexB This result is valid even when the trapping po- 
tential Uo is finite, which we demonstrate here. Since 
Eqs. (38) and (39) are valid in the presence of a fi- 
nite trapping potential, the transverse correlation func- 
tion can be expressed by wave functions in exactly the 
same form as that of vanishing trapping potential, up to 
Eq. (54). The wave functions, particularly those for core 
states, may be strongly affected by the trapping poten- 
tial, and may even become ill defined. An example may 
be the trapping of a vortex by a physical wire. Now, one 
may perform the same calculation of turning the area 
integration into line integrations, as done in Eq. (55). 
Since the trapping potential will not affect the superfiuid 
number density far away from the vortex, and since the 
circulation current is still zero at the vortex position, one 
then gets the same result as Eq. (55) in the presence of 
a trapping potential. 

The validity of Eq. (55) in the presence of a finite trap- 
ping potential implies that the transverse force is inde- 
pendent of the trapping potential Uq at the vortex center. 
In Eq. (29), the main function of the trapping potential 
is to specify the vortex position, a symmetry breaking in 
an otherwise homogeneous system. This is similar to the 
symmetry breaking by an infinitesimal field near a con- 
tinuous phase transition in statistical mechanics. Hence, 
it can be effectively taken to be zero, as we have explicitly 
done in the present paper. 

Next, we turn to the calculation of the superfluid num- 
ber density ps- At zero temperature, it is straightfor- 
ward. It is equal to the total fluid number density 
Pa = J2a.E^>o\^ai\x - xq\ -> oo)^, the number of 
Cooper pairs per unit area. At finite temperatures, there 
is a reduction of superfluid number density due to the 
backfiow carried by quasiparticle excitations. In principle 
one may directly calculate the current density together 
with the gap equation, or self-consistent equation, to find 
out Ps- This would be prohibitively difficult. Instead, one 
may proceed in the following manner: Far away from the 
vortex core, the current varies slowly. One may take the 
current to be locally uniform. Following the same way aa 
that in superfluid IIe3 using the backfiow contribution,LiJ 
the superfluid number density can be found as 



Ps{T)=po{l-Yo{T)), 
with the Yosida function Yq defined as 



(56) 



Yo{T) 



de 



,^£2+A^/feBT 



(eV^^^^^A'-^ + l)' 



which accounts for the quasiparticle excitations contri- 
butions. At the superconducting transition temperature 
Tc, Aoo = and Yq{Tc) = 1, the superfiuid number den- 
sity is zero as expected. This expression is the same as 
that obtained from the Londan penetration depth for a 
clean type-II superconductorB 

Using Eq. (55), the transverse term in the effective 
action, Eq. (34), is 



-- / dTdT'F^ir - r')(fa„(T) X Sxy{T'j) ■ z 



= i drdr' B {6xy{T) x (5i;i,(r)) • z 
= —i2TThps / dr 6xy{T) ■ At 



(57) 



with 



At = -ziSxy X z), 

which has the same form of the action for a charged par- 
ticle in a uniform magnetic field. The geometric phase 
or the Berry phase for the vortex moving along a closed 
trajectory F is 

Q ~ 2Trfips / dr 6iv{T) ■ At 

= 2n'hps (p d{Sxv) ■ At 

= -2TThpsS{T) 

with 5(F) being the area enclosed by F. The total trans- 
verse force on a vortex is then 

F = —2'Khps6xy X z . 

In view of the foregoing discussions, we may rewrite 
our general formulation, Eq. (34), in a more suggestive 
form. The effective action for the vortex is 



S. 



eff 



/ dri — i 2nhps5xv{T) ■ At 
+ - / dT'F||(T-r')|fe„(T)-fe4r')n (58) 



with At = ^{Sxy X z). The rewriting of Eq. (46) for the 
correlation function is 



^11 (^) 



1 



cosh 



dwJ{ijj)- 



(f-M 



sinh 



np 



and the rewriting of Eq. (47) for the spectral function is 
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•^(^) = I E -^(^^ - 1^" - ^«'l)l/a' - /J 



d32;<(a:)V,„Ho*a'(x) 



The thermodynamic hmit must be taken first to have a 
smooth spectral function, which is crucial for obtaining 
a finite vortex friction. 



IV. VORTEX DYNAMICS IN CLEAN LIMIT 

The example on extremely clean limit of fermionic su- 
perfluids is the superfluid He 3: The impurity concen- 
tration can be made to be smaller than 1 in 10^^. For 
superconductors, the impurity effect can, in principle, be 
made arbitrarily small, but no clear experimental realiza- 
tion has been reported yet. In view of this experimental 
situation, the discussions in this section are more relevant 
to He 3. However, from the methodological point of view, 
it is instructive to see how the formulation developed in 
Sec. HI works for such a clean situation. 



A. Extended States Contribution to Vortex Friction 

In this subsection we first calculate the extended state, 
quasiparticle and hole excitations, contributions to the 
vortex friction to illustrate the usefulness of the present 
longitudinal response formula. The formula, Eq. (47) 
or (49), is formally exact. However, for a given problem 
it is difficult to obtain an exact det ailed expression for 
friction, except in some rare cases.EjcJO Hence a WKB- 
type approximation will be used below. The responses 
of fermions, or electrons, governing by Hamiltonian dy- 
namics generates a finite friction for the vortex. 

At finite temperatures the extended states above (be- 
low) the Fermi level (the quasiparticles (holes)) are par- 
tially occupied. The vortex motion causes transitions 
between these states, which gives rise to vortex friction. 
The transitions between different single quasiparticle lev- 
els (^alVa^oHolvl/a') are considered here since they dom- 
inate the low-energy process. The quasiparticles are de- 
scribed by the eigenstates, Uq, and Va, of the Bogoliubov- 
de Gennes equation. Their behavior in the presence of a 
vortex has been well studied in Ref. pQ. We may take 



*a 



_ f Ua{x) \ _ 
Va{x) 



r^ikzZ 'ifi6-\-i{TzO /2 



Vl 



12-K 



-fir) 



(59) 



with r measured from the vortex position, is the az- 
imuthal angle around the vortex, L is the thickness of the 
superconductor film (the length of the vortex line) , and 
fo is the coherence length. In order to obtain a concrete 
form for the transition elements, we use a WKB-type 
solution for /(r) 



fir) 



1 

71 



^ Vg^-|A(r)r 



-*- + E 



n 



J^^.{k±{E)r) 



(60) 



with kp 



Here k±(E) = y'fc2 ± 2m^E^ - |A(r)|2/?^ 

fc? — k'^. The negative energy wave functions deter- 
mined by the Bogoliubov-de Gennes equation may be 
constructed according to Eq. (37) from the positive en- 
ergy ones. We will use the approximation that fcp ~ kp 
for the prefactor by assuming that the significant con- 
tributions come from the region near the Fermi surface. 
This WKB-type solution is valid when r is outside the 
classical turning point rt — \^\/kp. Here rt is the impact 
parameter. A WKB-type solution also exists inside the 
turning point. However, because it approaches zero as 
(rfcp)l^l/|/i|!, the contribution to the transition elements 
from this region is small, and will be set to zero. The 
transition elements are then given by 



d^xiu*^, (x) ( V^„ A)va, (x) + <, (x) ( V^„ A*)u^ix)) 







(61) 



Here Aoo is the value of |A(r)| far away from the vortex 
core. Physically, it means that if the classical quasiparti- 
cle trajectory is far away from the vortex core, it will not 
contribute to the vortex friction. The summation over 
states in Eq. (47) or Eq. (49) is replaced by 



E- E 



dEdE'- 



E 



E' 



2m 

^E^ - A2^ ^E'^ - A2j [l:' 

(62) 



after considering the density of states. 

Substituting Eqs. (61) and (62) into Eq. (49), 
and using the quasiparticle distribution function /„ = 
l/{e^^" + 1), the coefficient of friction is given by 






dE 



E^ 



1 



E^ - A^ cosh^ {pE/2) 



(63) 



The integral in Eq. (63) diverges logarithmically. It 
implies that the spectral function corresponding to 
the vortex-quasiparticle coupling is not strictly Ohmic 
but has an extra frequency factor proportional to 
ln(Aoo/?iw)- When tko is not very small comparing to 
Aoo, which may be realized when close to T^, we can ig- 
nore the logarithmic divergence in Eq. (63) by using the 
density of states for normal electrons to obtain a finite 
friction, i.e., replacing E'^/{E'^ — A^) with 1 in Eq. (63). 
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Close to Tc, the vortex friction approaches zero the same 
way as A^ , which is proportional to the superfluid num- 
ber density ps- When — ln(?ia;/Aoo) is large, we need to 
use Eq. (47) instead. Straightforward evaluation shows 
that in such a case 



■q 



Lm^CoAl^P 



1 



iTT^n^ cosh^ (/3Aoo/2) 



ln(Aoo/?itJc) . (64) 



Here lUc is the low-frequency cutoff. It is determined by 
the size of the system for a single vortex, and by the 
intervortex distance for a vortex array. 

We discuss briefly here the connection of our results 
to previous ones. The partial wave analysis has been 
performed for cuiasiparticle scattering off a vortex in a 
superconductor .E3 Though the phase shifts were obtained 
approximately, it is clear from the analysis that they are 
not all zero. Using the forraaJ |*ekttionship between the 
phase shift and the friction |H3i2aS the extended states 
have a contribution to the vortex friction, in accordance 
with our results. 

It should be emphasized that the logarithmic diver- 
gence conies from the interplay between the divergence in 
the density of states and the off-diagonal potential scat- 
tering. We can consider a situation where we physically 
create a pinning center to trap the vortex and guide its 
motion. In such a case the vortex has a diagonal poten- 
tial. If the scattering is dominated by the diagonal po- 
tential, e.g., by the trapping potential C/qj an additional 



factor coming from \ua\ 



will remove this logarith- 



mic divergence. The friction on the physical trapping po- 
tential will be finite even above Tc without the vortex, as 
indicated in Refs. |27|-p9| and ^, though the total trans- 
verse force disappears because ps = 0. This again shows 
the sensitivity of the vortex friction to details. 

The vortex friction from extended states exists for both 
clean and dirty superconductors at finite temperatures. 
Close to the transition temperature, it scales linearly 
with the superfluid density, and is exponentially small 
when T « Aqq. For intermediate temperatures T ~ 
Aoo, using ,^0 ~ h^kplmAoo and A''(0) = rnkp j -K'^fi' , 
T] ~ LhN{0)/^'^^/kBT. When the impurity potential is 
nonzero, there is an additional contribution to the fric- 
tion, to be discussed in the next section. 

We mention here that there is another type of low- 
lying excitation, phonons, which may lead to an addi- 
tional contribution to vortex friction. This type of exci- 
tation can be described by the phase dynamics of the gap 
function A, and has been ignored here by the assumption 
of an adiabatic following up of the gap function to the 
vortex coordinate. Based on general considerations w© 
expect that the phonon contribution is super-Ohmic.til 
Hence, it is asymptotically weaker than the (sub)Ohmic 
damping contribution from quasiparticle excitations dis- 
cussed above and the core state contribution to be dis- 
cussed below. 

The nonzero friction contribution from extended states 
found above is in accordance with the linear response the- 



ory in nonequilibrium statistical mechanics, where trans- 
port coefficients are related to the fluctuations near the 
equilibrium by the fluctuation-dissipation theorem. The 
fluctuations are completely determined by the underlying 
Hamiltonian. Here they are quasiparticles determined by 
the Bogoliubov-de Gennes equation. 

It should be pointed that the expression leadingto vor- 
tex friction, Eq. (47) or (49), are absent in Refs. 0,^ and 
O. To obtain a finite friction in those work, a finite re- 
laxation time needs to be inserted into-the denominator 
of the force-force correlation function,E2l or into the de- 
nominator of theJNambu-Gor'kov Green's function at a 
convenient point p which at the same time leads to the 
reduction of the total transverse force. As discussed in 
Sec. II. B, such a procedure should be avoided. 



B. Core v.s. Extended States Transitions for Total 
Transverse Force 

There are various ways to express the transverse cor- 
relation function in Eq. (36) or (52), with emphasis on 
different aspects of the transition elements. In Sec. III.C 
we have shown how to obtain the total transverse force 
from the consideration of extended states. In this sub- 
section we show that it can also be obtained from the 
consideration of core states, even some combination of 
both types of states. 

Explicitly, we may evaluate the transverse response di- 
rectly from Eq. (36) or (52). We will show that the total 
transverse force can be expressed as contributions from 
only core to core transitions, or from core to extended 
states transitions. For a clean superconductor we replace 
VxQ -^ —V. We define following symbols for the transi- 
tion elements; 

a{l)a'{l') = -ih{fc, - fa')^z ■ f d^xcfx' 

*t^(a;')V*a'(x') X V*^,(x)*„(x) , (65) 

which groups the transition elements into core state to 
core state, core state to extended states, extended state 
to core states, and extended state to extended state tran- 
sitions. Here I = c,e represents the core c or extended e 
states, and a represents other indices: kz, fJ-- For exam- 
ple, a{c)a'{c) represents the elements in Eq. (52) when 
^I'q and ^Q-' are both core states. More explicit exam- 
ples will be given below. ^From Eq. (52) and (65), the 
summing over these transition elements as well as over a 
and a' gives B, 

B = ^ [a{c)a'{c) + a{c)a'{e) + a{e)a' (c) + a(e)a'(e)] . 

Q.a' 

(66) 

First, we note that for a core state, the sum of its tran- 
sition elements to all other states is zero: 
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Y^ [a{c)a'{c) + a{c)a'{e)] = . 



(67) 



In fact, we have already obtained Eq. (67) and used 
this identity earher in Sec. III.C from Eqs. (52)-(55) 
to exclude the core state contribution to the circulating 
current far away from the vortex core in Eq. (55). In 
Eq. (52) both extended and core state contributions are 
there. Then we summed over all the ^cC to reach Eq. 
(53) or (54), because they form a complete set. Thus in 
the last expression of Eq. (52) inside the sum over a, if 
\E'q is an extended state, its contribution to the transverse 
response is actually J2a' [a(e)a'(c) + a{e)a'{e)], because 
all the a's have been summed over. The same procedure 
applies if ^a is a core state. In Eq. (55), we have shown 
that the area integral can be converted into a line integral 
and we can choose the loop large enough compared to the 
core size. If '^a is a core state, its contribution to Eq. 
(55) is zero. Thus we conclude the validity of Eq. (67). 

With the aid of Eq. (67), the transverse correlation 
can be expressed in the following two additional forms: 



B = Y^ [a{e)a'{c) + a{e)a'{e)] 

a, a' 

= J2[a{e)a'ie)-a{c)a'ic)] 



(68) 
(69) 



In reaching Eq. (69), we have used the identity 
a(c)a'(e) — a{e)a'{c). Here Eq. (68) can be reduced to 
Eq. (55) after using the completeness of the eigenfunc- 
tions, as discussed above. Next, we present more detailed 
discussions on the core-core and extended-extended tran- 
sition element contributions. 

To be specific, we consider the zero-tcmpcrature case. 
For the core states, because of the topology the energy is 
uniquely determined by /i. The only transition elements 
contributed to the transverse correlation function F±_ are 
between states /i = i^: 



J2(^ic)a'{c) 



a, a' 



~i2h 



(/- 



,=-. f,,=,^)-z- I d'x I d'x' 



hl,ix')V,,M>.ix') X V,„vl,|(x)*_i(. 

\ 2 ^ 2 ^ 



(70) 



The additional factor 2 accounts for the transition from 
the /i = 1/2 state to the /i' — —1/2 state, which gives 
the identical contribution. The transition element in Eq. 
(70) may be expressed as 



I^'^I 



Sx / d?x'Ul,{x')^.,,^i{x') X V,„*t(x)*_i(x) 
{a_i 1 +6_i iVa* 1 1 +6* 1 i) 

V 2^2 2'2y V ^2-2 "2^2/ 



-izSk^.k' 



^{E-i-E^' 



d-i 1 = — 

2 ' 2 



rdr\At tdr)fUr) 



+t^dr)UAr) 



and 



h_L 1 

2 ' 2 



dr\A\ fl^_dr)f-i(r) 



-nAr)f+Ar) 



which follow the definitions in Eqs. (A10)-(A15) in the 
Appendix. 

Now we evaluate Eq. (70) explicitly. For the deep core 
states in clean superconductors E = — /ieo with eq th^ 
core level spacing, fu=-i. = and /„^i = 1. The rela- 
tion between the energy of core states and the quantum 
number jj, in our case is different from the one in Ref. |30| 
in sign because we are considering a vortex with positive 
vorticity. The wave functions for deep core states take 
the form 



H(^) 




e'(^+*'^Vi(fcFr) 



-r/io 



(71) 



This leads to a_i/2,i/2(S, E') « and 6-i/2,i/2(£^, E') 

Aoo/^. 

Using Eq. (71), it is straightforward to show 
5 • E(/a. - f,')jd\jd\' {¥^{x')^^,.{x')x 
V*;,,(a;)vl'^(x)) ^ikj 



k\. 



Therefore, 



^a(c)a'(c) = ^ = B 



2 



(72) 



where 



The last equality is due to the fact that in 2D, the elec- 
tron density n,, — k'^p/2'K. The additional factor 1/2 ac- 
counts for the pairing. The conclusion is that at zero 
temperature the sum of the core-core state transitions 
alone gives rise to the total transverse force. It corre- 
sponds to the fact that the core-core state transitions are 
a local and differential form of the geometric phase, and 
the Berry phase is the global and integral form. In the 
next section we will show that Eq. (72) is unchanged in 
the presence of impurities. 

Here we wish to point out an interesting feature ex- 
plicitly manifested in Eq. (72): The transverse response 
is insensitive to the size of the system, because the core 
states are exponentially localized. This implies that the 
thermodynamic limit is not important for the total trans- 
verse force. We attribute this feature to the topological 
constraint on the transverse response, corresponding to 
the well-known fact that the Berry phase exists for a dis- 
crete energy spectrum. 
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There are two more interesting results which follow Eq. 
(72). At zero temperature, using Eqs. (69) and (72), we 
have 



a(e)a'(e) = 2B 



(73) 



This implies that for a fermionic superfluid, the sum of all 
extended state transitions lead to twice that of the total 
transverse force. Combining Eq. (67) with Eq. (72), we 
have 



B = -Y,a{e)a'{c) 



(74) 



which shows that the core-extended state transitions can 
also be used to calculate the total transverse force. We 
believe that this property has been explored before in 
the case considering the contributions from states whose 
energies are around Aoo, the interface between the core 
and extended states. lJ 

In the literature, after a transverse response equiv- 
alent to Eq. (36) or (52) was reached, it had always 
been assumpd_that only one core-to-core state transition 
contributesJjH However, as we have found out, core to 
extended state transitions are of the same order. The 
above discussions show that there are many equivalent 
ways to compute the total transverse force. Because of 
the topological constraint, the total transverse force can 
even be evaluated by partial summations of the transi- 
tion elements, expressed by Eqs. (68), (69), (72), (73), 
and (74). This is completely different from the computa- 
tion of the longitudinal force (the friction) where a partial 
summation contributes only a part of the total friction. 
The demonstration in this subsection suggests that the 
alleged cancellation between the core spectrum flow con- 
tribution and the Berry phase counting is a consequence 
of the combination of double counting, treating core and 
extended contributions to the transverse force as different 
quantities, and the misuse of the relaxation time approx- 
imation in the force- force correlation functions. 

To briefly summarize this subsection, we have formu- 
lated the transverse response in terms of transitions be- 
tween core state or between core and extended states. 
Equivalently, we have also formulated it in terms of a 
summation over the extended state contributions. In a 
clean and neutral superfluid at finite temperatures, the 
total transverse force is given by the product of the su- 
perfluid number density, the Planck constant h, and the 
vortex velocity, though the vortex friction exists. 



V. EFFECTS OF IMPURITIES 

A. No Effect on Total Transverse Force 

The presence of impurities is unavoidable in supercon- 
ductors. In this subsection we consider this realistic sit- 
uation of the influence of the impurity potential to the 



transverse force on the moving vortex. In Sec. IV. B we 
have shown that the transverse correlation function can 
be evaluated by either considering the extended states or 
by considering only the core states in a clean supercon- 
ductor. The same also holds in the case with impurities, 
as we will demonstrate below. We first give a formal 
demonstration from the counting of individual state con- 
tributions, then explicitly consider the core state tran- 
sitions, to pave the way for the core state contribution 
to vortex friction. The robust conclusion is that random 
impurities do not affect the total transverse force. 

It is more convenient to change the gradient from V^,;, 
to V when there is an impurity potential V{x) involved. 
Applying Eqs. (42) and (43) to Eq. (52), the transverse 
correlation becomes 

h ■ (*i(.T')V*„-(a:') X V^l{x)^^{x) 



d?x / (fx' {^l{x')Sin'{x')^o.,{x')y. 
^i,{x)V7i'{x)-^o.{x)] . (75) 



Other terms are identically equal to zero after summing 
over a, a'. 

We first show that the second term in Eq. (75) is zero 
after the impurity average. To be concrete, we expand 
the wave function $„ in terms of eigenfunctions of Tio, 



*a = X! ^"^ 



(76) 



Here Xa-y — aa-ye'"'^'''' , and Oq,^ and ipa-f are the modulus 
and phase of the expansion coefficients. We remind the 
reader that Hq has dependence on the impurity potential 
because it includes the smooth part of the self-consistent 
potential A. Because of the normalization requirement, 
the coefficients {aajC^''^"'^} form a unitary matrix, 

a 

The second term of Eq. (75) is now 

d^x I (f-x'-^l{x')\in\x')-^a,{x') X 
^l,{x)Vn'(x)^c.(x) 

^Q7'^a'7'^a'7j ^071 ^ i i X 

7:7':7l:7l 



/ , ""07 

7 



d?x / (fx'(j)\{x')V-H'{x')(l)Y{x') X 



4>\,y)vn'{x)<i>^,{x). 



(77) 
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^From the random-matrix theoryJ2a the phase {1^90:7} are 
random numbers. The impurity average makes 7 = 71 
and 7' — ^'y- Under the assumption, i.e., the core size 
much larger than the average distance between impuri- 
ties, the averaging over impurities restores the homogene- 
ity and isotropy of the spatial space. This implies that 
all the odd power of V7i' will be averaged to zero. Since 
each term in Eq. (77), 



d^x I d^x'<p\{x')\'H'ix')(l)Y{x') X (/)t,(a;)VH'(a;)0^(.T) 
z I d^x I d^x' U(f>\{x')'^n'{x')(l)Y{x'))^x 



( 



^i,{x)^n'{x)<f>^{x)) - {cf>!^{x')^n'{x')<f>^,{x')) X 

(^^l,{x)Wn'{x)cl>^{x))^} , 



consists of the elements of odd power of x and y compo- 
nents of V7i', the second term of Eq. (75) is zero after 
the impurity average. 

The average transverse correlation is then, following 
the same procedure from Eqs. (52)-(55), 



B = (j) dl • j 



(78) 



Here the average current 



ih 



j = - < y ^ [/a<VM„ + (1 - /a)WaV<] > +C. 



2 



< "17 > {[faU*'Vu^ 



+ {l-U)v-yVv*]+c.c.} , 

and < ... > stands for the impurity average over the 
expansion coefficients. In the limit r — > 00, we have 

j = -yl]<aa,7> {[/"l"7l' 

a, 7 

+(i-/„)|«7|2] +c.c.}ve 

At zero temperature, we have 



J2<alj> KQx - Xo\ -^ oo)|^ 



a,Eo.>0 7 

= ^ < \va{\x - xa\ -^ oo)\'^ > 

q,-Eq>0 

= Po • 

The above second equality is the Anderson theorem, 
in that nonmagnetic impurities do not affect the den- 
sity of states near the Fermi surface, hence there is no 
effect on the superconducting transition temperature. 
We will come back to this point after the discussion 
of the impurity effect on vortex friction. This result 
may also be reached from the envelope wave function 



argumentationBa Therefore Eq. (55) also holds in the 
presence of impurities. Thus we have shown in general 
that the transverse correlation is not influenced by impu- 
rity potentials. Physical understanding of this result is 
straightforward: There is no average circulation current 
associated with impurity potentials. 

Next, the insensitivity of the total transverse force 
to random impurities will be illustrated by a different 
demonstration. We evaluate core state transitions with 
impurity potentials, a part of the first term in the right- 
hand side of Eq. (75). First, we will explicitly consider 
the total transverse force from core states contributions 
with a weak impurity potential. Because of the factor 
(/m ^ ffj.') and the selection rule in the transition ele- 
ments, deep core states are the ones making the main 
contributions to the total transverse force for temper- 
atures well below Tc- When the impurities are weak 
enough, for deep core states ^^t — J2iy ^tiv^^'^'"' 4>v the 
expansion coefficient a^^ is large only for the neighbor- 
ing states around /x, which are also deep core states. In a 
clean superconductor, core states are uniquely specified 
by the azimuthal number u. Therefore, with weak im- 
purity potentials, we may only consider deep core states 
close to the Fermi surface, and ignore the mixing of deep 
core-level states with the extended states in Eq. (75). 

Substituting Eq. (76) into Eq. (75) and using Eq. (71) 
to calculate transition elements among core states {4>^}, 
we have 



d^x / d^x''^Ux')V'^^,,{x') X V*)^,(a;)*^(x) 



-i-Pi^v+iV^'u'-iVf^i^i +i'P^ui 



I'm' Ml m{ 

Here we have used 



d^x(j)l{x)V(j)^> (x) « itc {±x + iy) 5^' 



(79) 



v±l 



with tc = kp jl. If /i and /i' are interchanged, there is a 
sign change in the right-hand side of Eq. (77) . Including 
the factor (f^ — ffj,') and summing over /i and fi' , core-core 
state transitions become 



E 

z ■ 



iU-U') I d^x I d^x' 



■[^l{x')V^^,{x')xV^l,{x)^^{x) 



^-.^.„+i^,,.,-..,,„,+.^.„, 2(T)*t?<J.',.±i5.;,.,±i. (80) 

Because of the randomness in the phase factor of 
e ^ '1 , tlie dommant contribution 

comes from states v — vi^v' — v[. Equation (80) be- 
comes 
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A',M' 



(ji-i jfj.') 7 ^ 



2 2 



,2(T)it%, 



v±l ■ 



For a given distribution of Cq,^ this summation can 
be evaluated. For the purpose of demonstration, let us 
assume a simple distribution centered at a: a^^, = 1/(2/+ 
1) when l?^ — /i| < I, a^i, = otherwise, and consider only 
zero temperature. Here 1 << / << A^o with Nc being 
the total number of core states. With this assumption, 
each of the original states is spread into a band of 2/ + 1 
states around it when impurities are present. Equation 
(79) takes the value 



n,n' 



[Jfi jfj,' ) 



E 



|f— /j,|</,|i/'— /j,'|</ 



{21 + 1): 



:2{T)it%',^±i 



We note that for the pair of states /i = Tl/2, /i' = ±1/2 
closest to the Fermi surface their contribution is reduced 
by a factor of 1/(2/ + 1)^. However, all the states within 
the energy shell \E^\ < ho near the Fermi surface con- 
tribute now. In order to have a nonzero contribution at 
zero temperature, we have E^Ef^i < 0, one above and 
one below the Fermi surface. Because of the restriction 
of the band distribution and the selection for v' and v, 
we have an additional constraint on the E^j, and Ep,/ : 
\Efi\, \Ef^i\ < I. The net contribution is 
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which is approximately the same value for the clean su- 
perconductor. The factor of 1/2 is due to the approxi- 
mation for tc using same value for all transitions between 
\x and /i ± 1. One can check that the above result also 
holds for Gaussian distribution of Cq^. This completes 
the discussion of the first term in Eq. (75) in the weak 
impurity potential limit. 

Although at zero temperature in a clean superconduc- 
tor only the core states closest to the Fermi level con- 
tribute to Eq. (52) or (66), the transition elements of 
other states are not small. Their contributions cancel 
each other completely. With impurities, more states than 
those closest to the Fermi level give contributions to the 
transverse response. These contributions from other core 
states restore the transverse response to its original value 
of a clean superconductor. In the calculations with the 
relaxation time approximation, the reduction of the con- 
tribution from the two stateSr-dosest to the Fermi level 
has been taken into account .Q't3 The contribution due to 
other core states, which arises after introducing impuri- 
ties, has not been included in those calculations. 

Next we consider the dirty limit, and we will again, 
make use of arguments in the random matrix theory.LJ 
We assume that there is no mixing between the core and 
extended states. In the weak impurity potential limit, it 
is not difficult to justify this assumption: The band width 



in Eq. (81) caused by impurities is much smaller than 
A. In the dirty limit, the number of core states remains 
the same, since the energy gap away from the vortex core 
remains the same. Hence, there is a conservation of the 
number of core states, because of the topological nature 
of the vortex. We also note that there is no degeneracy 
for the core states, in contrast to the extended states. In 
additional, impurities do not cause an additional viola- 
tion of time-reversal symmetry. For those reasons we do 
not expect that they would mix two topologically distinct 
types of solutions, the core and extended states, of the 
Bogoliubov-de Gennes equation. 

With increasing impurity potential strength, eventu- 
ally any core states VP^ in Eq. (79) consist of all the 
core states of Tio: ^/j = ^i, Xni' 4'u- Here the summation 
^j^ runs over core states only. The total number of core 
states does not change after introducing impurity poten- 
tials, because the core-level spacing for (l)i, only depends 
on the values of Ep and Aqo, the value of A far from the 
core. A specific approximate realization may still be in 
the form of the band distribution as given in Eq. (81), 
with I r^ Nc, the total number of core states: 
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It is widely spread for deep core states |/i| << Nc/2. One 
may check that Y,„ IXm^^P = 1 and M = Li. IXm^'P'^- The 
latter corresponds to the requirement for the energy spec- 
trum < E^, >^< *^|Ho|*M >= E^ \Xt^,^?E^. The con- 
2) dition of ^ IX/^i^P = 1 is only approximately satisfied: 
we have found that E^ IXm^P ~ \ ^^{^l)l{{^cl2f-v'^) 
which gives In 2 when v << Nc/2 and InA^c when 
1/ ^ Nc/2. 

In this limit, Eq. (79) becomes 
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Its average value is ikp/2 at zero temperature, approxi- 
mately the same as in clean superconductor, by using the 
distribution function given in Eq. (82). 

We have gone into great detail to calculate the total 
transverse force from core state transitions in the pres- 
ence of impurities. Indeed, there is a reduction in tran- 
sition amplitude between any pair of neighboring states. 
Nevertheless, the summation over all possible core state 
transitions restores the total transverse force to its value 
in the clean limit. Hence, the impurities have a neg- 
ligibly small effect on the total transverse force from 
both the core state transitions consideration and the ex- 
tended state counting, though the friction contributions 
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are strongly affected by impurities, to be discussed in the 
next subsection. 

In addition, we check the self-consistent condition with 
respect to A here and show that they are satisfied for 
our choice of Xp.i'- Because in Tin we have already as- 
sumed that the profile of A is the smooth part of the self- 
consistent potential with impurity potentials included, 
we need to make sure that the decomposition of the 
eigenfunction does not introduce an extra term to the 
self-consistent potential. We have 
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where u^,v^ are the components of core eigenf unctions 
{01/}, and (e) denotes the extended states. In the last 
equation, the summation is split into those of core states 
and extended states. Using a distribution function p for 
X such as defined by Eq. (82), the average core state 
contribution to A is the same as the one calculated by 
using {0^}, 

The extended states also need to be self-consistent. We 
assume the impurity strength is strong enough to mix 
the core state on the scale of 1/^ but too weak to cause 
extended states distortion on the scale of l/kp. Then 
for extended states the distribution of Xa-y is a function 
p{Ea — E^). It is straightforward to check that the ex- 
tended state contribution to A is the same as that of 
those calculated by using {(fiu}- 



B. Impurity Contribution to Vortex Friction 

In the present of impurity potentials, there are two 
kinds of contributions to the friction. The extended state 
contribution remains basically the same as what we have 
discussed before. The main difference is that the co- 
herence length in Eq. (61) will change when impurities 
are present. We will give a brief discussion here. Since 
the density of states remains unchanged, we only need 
to evaluate Eq. (61) again. The transition elements are 
given by 
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Here x is the expansion coefficient in Eq. (76), ^ is the 
coherence length in the presence of impurities, I is the 
angular index of the state, and E is the energy of the 
state. With a given distribution of x such that the ex- 
pansion coefficient is confined to the neighborhood of its 
original energy, it can be shown that the extended state 
contributions to the friction remain unchanged, except 
the change of the coherence length of ^o — > C- 

In the presence of impurities, the core state energy lev- 
els are no longer monotonically arranged according to az- 
imuthal number or the angular momentum. In addition, 
it may become quasicontinuous under the impurity av- 
erage. The mixing caused by impurity potentials makes 
it possible to have transitions into energetically nonadja- 
cent core states, as discussed in the previous subsection, 
as well as into energetically nearly degenerate core states. 
Thus the core states can give another contribution to 
the vortex friction, similar to the residual resistance in a 
metal. 

First, we consider the weak impurity potential limit. 
We assume the effect of impurities is not so strong such 
that we can treat their influence on core states perturba- 
tively. Using Eqs. (42) and (43) we obtain 
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For the leading-order contribution, we only need to use 
the unperturbed \E'^ and A. The first term in Eq. (86) 
will not give any contribution to the dissipation because 
of the discreteness of {E^> — E^j) and the factor 5{fito — 
|£'^t — £'^'1) in J{u!). After summing over ^ and /j,', this 
term will become terms of 5{hu> ± eo), which will not 
give any dissipation. The contribution to the dissipation 
comes from the second term. 

We may assume the impurity potential has a length 
scale small compared with the coherence length so that 
we can describe it by a delta potential V{x) = ^^ Vq5{x— 
Xi). We have 
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with the impurity concentration rii. Under this assump- 
tion we will make connections to the normal-state trans- 
port parameters. Note that for normal states the elec- 
tronic transport relaxation time and the electmii scatter- 
ing cross section have the following relationsiEZl 

tZ} = n.,VF(Jtr 
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with 



Here 
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with q = k — k', is the angle between k and k', and 
vp — hkp/m. For our choice of impurity potential, r^^ 
can be calculated, 
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We emphasize here that the electronic transport relax- 
ation time Ttr is directly determined by the Hamiltonian 
ofEq. (22). 

Expressed in tj,., the spectral function now becomes 
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Here eo = A^^/Ep is the core-level spacing, and Ep ~ 
mVp/2. In Eq. (87) kpL/i-Kto is the approximate density 
of core states near the Fermi surface. It appears with the 
factor 2fiLd because 

Y,5{huj~\E^~E^.\)\f^-f^,\ 

= j dE^dE^,5{riLo - \E^ - E^,\)\f^ -U'\y- nl{E) 

= 2hujnl{E) 

with nc{E) Ri kpL/{2'!Teo). The scattering time Ttr is 
linked to the residual resistivity by 
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and can be measured independently. 

The above spectral function J{uj) gives the vortex fric- 
tion in the weak impurity limit 



V 



3 mrie^Tr^'^ L) 



Ttr 



It has a simple interpretation. For a normal electron 
moving in the metal, the friction is simply m/Ttr- Equa- 
tion (88) can be interpreted as that in the weak impurity 
limit, the friction for a vortex is the friction for each elec- 
tron times the total number of electrons inside the core. 

This vortex friction increases with impurity concen- 
tration and strength. We will show that this increase 
eventually saturates in the dirty limit. Using Eq. (76), 
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eigenfunctions of Tig- 
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With the distribution function given in Eq. (82), the 
average value 
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Here the total number of core levels is 
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Finally, the spectral function is 
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which gives the friction per unit length 
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This result is similar to what is obtained in Ref. |l|. Hence 
its microscopic base has been provided. In the low- 
temperature limit, the magnitude of the vortex friction 
is smaller than the total transverse force by a factor of 
A^/Ep. 

In the above derivation, we have ignored the local- 
ization effect which suppresses the density of state, or 
the superfluid number density. We justify our assump- 
tion here. There are three energy scales involved in the 
derivation of vortex dynamics, the Fermi energy Ep, 
the energy gap Aqo, and the core level spacing A^/Ep. 
The effect of impurities on vortex, dynamics is believed 
appreciable at TtrA^/fiEp < Ip and the equality of 
Eqs. (90) and (87) suggests that the impurity starts 
to be effective at Ttr Aoo /h {Aoo /EpY ^ 1- They in- 
dicate that the impurity effect on vortex friction occurs 
at a rather weak level, determined by the smallest en- 
ergy scale in the problem. The dirty limit is given by 
Aao/Ep < TtrAao/h < 1. The localization effect is only 
pronounced in the extremely_d,irty limit, the localization 
regime, when TtrEp/fi < 1.E3 Because Aao/Ep « 1, 
away from the localization regime the suppression of den- 
sity is indeed negligible. The unsuppressed electronic 
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density applies, and the present results are valid well into 
the dirty limit of the superconductors. 

To summarize this section, we have shown that the 
total transverse force is insensitive to impurities by two 
different methods, but the additional core contribution 
to the vortex friction arises. For a weak enough impurity 
potential, a perturbative calculation leads to the core fric- 
tion proportional to the normal-state resistivity. In the 
dirty limit the core friction contribution saturates to a 
value determined only by the energy gap and the Fermi 
energy. 



VI. COUPLING TO ELECTROMAGNETIC 
FIELD 

Now let us discuss a superconductor when the pen- 
etration depth is finite but still much larger than the 
coherence length. The Lagrangian is given by 



1 fh 



Lbcs = ^ -(Plix, t) I hdr - MF - eAo + -— I - V 
--^A^ + V{x)\^,{x,t) 



-gipl^x, t)^\ {x, t)V;| (x, r)V'T {x, r) 

+ ^(i;2-fB2)-feAono, (91) 

here —chq is the charge density of the ionic background. 
The coupling to the electromagnetic field is in the usual 
minimum coupling form. The fermionic degrees of free- 
dom can be integrated out to give 
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with 
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Here 



Zm \ I c 



H* = -eA„ + l-f^y + -^a\ - mf + V{x) , 

and c is the speed of light in a vacuum. Variation with 
respect to Aq and A gives 
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with E = -VAo + \dA/dT and B = V x A. Here 
en and ej are the electric charge and current densities. 
They should be obtained through the electronic Green's 
function. In the Lorentz gauge, the|-equations for Aq 
and A from the above equations arej^ adapted to the 
imaginary time here. 
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Assuming that for a static vortex at Xv the vector po- 
tential is A(x — Xy), then for a slow moving vortex, 
the correction 6A to the vector potential A from A = 
A(2; — Xy{T)) starts from second order in Xv(t) and can 
be ignored. This can be directly demonstrated from Eq. 
(96) to the leading order in SA: 
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The same is true for the scalar potential Aq. For our 
purpose of keeping to linear order in Sxy we may use 
A = A(a; — Xy{T)). 
Now we expand 

Al{x - Xy{T)) = (1 + SXy{T) ■ V^jo) Al{x - Xy{T)) , 

with / = (0, X, y, z). The effective action for the vortex is 
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with 
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where both E and B are substituted by the stationary 
values calculated from ^o = and A = A{x — Xy{T)). 
Equation (97) is in the same form as Eq. (36). There is 
then no change of the transverse response from the su- 
perfluid, reflecting the fact that the canonical momentum 
of the superfluid is not changed by the coupling to the 
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electromagnetic field. A similar conclusion has also been 
reached by others phenomenologicallyO 

However, there are relativistic corrections to the so- 
lutions (Eqs. (95) and (96) ) of the Maxwell equations 
due to the motion of the current and charge sources as- 
sociated with the vortex. They are determined according 
to the Lorentz transformation of the four vector formed 
by the scalar and vector potentials, or equivalently, the 
four vector by the charge density and current. Those rel- 
ativistic corrections give rise to additional terms in the 
action and can, in principle, contribute to the total trans- 
verse force in vortex dynamics. The relevant term in the 
effective action takes the form 

dr / (fix -Siv ■ A{x — Xy) [n{x — Xy) — tiq] , 
J c 

arising from the relativistic correctijOa to the scalar po- 
tential, the Aharonov-Casher phaseO This contribution 
is due to the interaction between the moving magnetic 
flux carried by a vortex to the electric charges of both 
conducting electrons and the background charges. Since 
the charge neutrality condition is maintained in a su- 
perconductor, this contribution is zero, as also been no- 
ticed in Ref. E2I Hence, we do not need to consider it 
here. Other relativistic corrections do not affect the to- 
tal transverse force, and the rest of the terms have the 
same structure as the uncharged superconductor with the 
replacement of Eq. (98) to Eq. (29). Therefore, all the 
steps from Eqs (24)-(36) remain unchanged. We arrive 
at the same expressions (Eqs. (35) and (36)). There 
are two differences, however. First, the Bogoliubov-de 
Gennes equation includes a vector potential, which can 
be served to generate a vortex, not by a rotation of super- 
fluid. Second, we have now a compelling physical reason 
to neglect the phonon (density fluctuation) mode com- 
pared to the neutral case, because it is the plasma mode 
with a big energy gap. 

Let us discuss the effect of including the vector poten- 
tial in the Bogoliubov-de Gennes equation to the final 
results. For the extreme type-H superconductor the vec- 
tor potential near the core is Ag = -^rhQ, here h^ is the 
magnetic field along the vortex line. When the pene- 
tration depth is large, /lo is small. For small r, when 
solving the vortex core structure, Tic/2er > Ag and we 
can safely ignore A0 = 2^^- The core structure is in- 
sensitive to the vector potential in the extreme type-II 
case. Because the total transverse force can be expressed 
in core state transitions, it is insensitive to coupling to 
electromagnetic field. Equivalently, when expressing the 
transverse response in terms of the summation over ex- 
tended states, Eq. (55), we need the large-r behavior 
of the Bogoliubov-de Gennes equation. When r >> A, 
Ag ^r Q^ so that the the coupling to the electromagnetic 
field will not influence the results in Eq. (55). In the pres- 
ence of impurities, the vortex friction is also insensitive 
to the cqjfipling to the electromagnetic field, as noticed 
long ago 



To summarize this section, in a charge neutral extreme 
type-H superconductor, the vortex dynamics is the same 
as that in an uncharged BCS superfluid. This is a known 
result, but we have sketched how to obtain it within the 
present formulation. 



VII. EXPERIMENTS 

A. Transport Measurement 

It had been assumed that the forces on a vortex could 
be extracted from transport measurements. Let us first 
review this apparently plausible proposal and discuss 
ideas which are crucial to the understanding of transport 
measurements. Considering our previous derivations, the 
equation of motion of the ith vortex takes the form of the 
Langevin equation similar to that of a charged particle 
in the presence of a magnetic field: 



= qvPsh (Vs - Yi 



rjYi 



Ef. 



(99) 



>f 



Here q„ = ±1 represent different vorticity. The total 
transverse force —qyPshr and viscosity rj are the ones 
we have calculated in the previous sections. In addition, 
there are a fluctuating force f related to the viscosity 
by the fluctuation-dissipation theorem, a pinning force 
Fpin, and we should also include the forces due to other 
vortices F^ because of vortex interaction. Here we have 
explicitly written out the external current term in Eq. 
(99), though in a real situation its effect is always through 
the rearrangement of vortices in the superconductor. The 
motion of vortices is a genuine many-body problem. A 
general exact solution does not exist. 

Equation (99) may be solved after a drastic simpli- 
fication by ignoring the pinnings. This is equivalent 
to the situation that a perfect vortex lattice is sliding 
through the sample. Together with with the Josephson 
relation, we can determine longitudinal and transverse 
resistivity for superconductors. The Hall angle, defined 
as 6Haii = t&n^^ {pxy / Pxx) , is nearly 90 degree for almost 
all situations. However, in transport measurement, most 
of the samples show a small Hall angle and some show a 
sign change in the Hall angle upon entering the supercon- 
ducting state. This simplified model certainly disagrees 
with experiments. 

Now let us consider whether or not this simplification 
can be made by considering the magnitude of pinning. 
The equation of motion of a vortex in a superconductor, 
Eq. (99), has the form of a particle with zero mass in a 
strong magnetic field. Because the kinetic energy is zero, 
it is always confined to a local energy minimum in space, 
formed by vortex interactions and pinning potentials. An 
applied current tilts the potential and the vortex moves 
by thermal activation. If the applied current is so large 
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that no minima due to pinning and vortex interaction ex- 
ist, then indeed we expect a large Hall angle. However, in 
order to have a truly free vortex flow, the current should 
be large enough to overcome the largest pinning poten- 
tial, the edge pinning. The current needed is on the order 
of 10^ A/cm^, which is too large to be relevant to exper- 
iments. Therefore, in the real experiments, the vortices 
must be helped by their many-body interactions to over- 
come this energy barrier. We need to consider transport 
measurement by solving the lattice structure formed by 
vortices and by what mechanism their transport is made 
possible. It has been quantitatively suggested that vortex 
many-body effects can be responsible for the Hall effect, 
in agreement with recent experimental indicationsJlj 



B. Direct Measurement of Total Transverse Force 

The total transverse force on a moving vortex in thin 
YBCOpQlms has been directly measured via a mechanical 
device.O This experiment used a small vibrating magnet 
mounted above the center of a superconducting film to 
generate moving vortices in the film. The vortices fol- 
low the motion of the magnet for samples with less twin 
boundaries. The experiments were performed on those 
samples. The force was measured by measuring the mo- 
tion of the superconducting film in response to the vi- 
bration of the magnet. The experimental results have 
provided a qualitative confirmation of the insensitivity 
of the total transverse force to impurities. 



C. Measurement of Friction 

In a rf resistance measurement, the vortices are mov- 
ing around their local minima, rather than over potential 
barriers in a DC measurement. In such a case, it is pos- 
sible to observe intrinsic friction after using a potential 
to describe the periodic vortex interaction and making 
some assumptions about the pinning. The rf resistance 
was analyzed early with a vortex dynamics model with- 
out transverse force.c3 In order to compare with our the- 
ory, the total transverse force needs to be included. We 
will not go into any detail other than to suggest this pos- 
sibility. 



VIII. CONCLUSIONS 

We summarize here what we have achieved in the 
present paper. With respect to the microscopic deriva- 
tion, we have developed an influence functional formula- 
tion started from the BCS theory This formulation has 
allowed us to discuss several difficult questions regard- 
ing vortex dynamics. One question has been whether 
the total transverse force originates from core states, ex- 
tended states, or from both. This question is unique to a 



fermionic superfluid because of the vortex core structure. 
We have shown that the total transverse force can be 
calculated equivalently by considering exclusively tran- 
sitions between core states, by transitions between core 
and extended states, or by counting contributions from 
extended states. The total core-state transition contri- 
bution to the total transverse force is shown not to be 
affected by impurities when calculated by using random 
matrices instead of the relaxation time approximation. 

On the thermodynamics and statistical mechanics 
level, we need to consider the increase of the superfluid 
kinetic energy associated with the increase of superfluid 
momentum due to the vortex motion. This kinetic en- 
ergy needs to be provided from somewhere. If there are 
no normal fluid and no impurities, this kinetic energy is 
provided from the work done by the external trapping 
potential on the vortex. When either the normal fluid 
or impurities, or both, are present, there is a question 
whether or not a vortex can extract the internal energy 
from the normal fluid or substrate, which carries entropy, 
and can transfer it into kinetic energy of the superfluid, 
which carries no entropy. If not allowed, the increase of 
superfluid kinetic energy due to vortex motion can only 
be provided by an external force and the transverse force 
on a vortex cannot be reduced by the normal fluid or 
random impurities. We have discussed this question and 
demonstrated that thermodynamics gives a powerful con- 
straint on phenomenological models of vortex dynamics: 
The total transverse force cannot be reduced. 

We have located the source for contradicting theoreti- 
cal results in the two pictures: the use of the relaxation 
time approximation in the force calculation. This prob- 
lem is rather subtle. It is well known that the relaxation 
time approximation has been used successfully in some 
applications, particularly in calculations of conductivity 
or mobility, where the average velocity is computed un- 
der a given driving force, through velocity-velocity cor- 
relations. However, the special feature of vortex dynam- 
ics is that it belongs to the same category as resistiv- 
ity or friction formulas in transport theory, where the 
average force is computed with a given velocity. The 
direct calculation of resistivity is known to be difficult. 
To derive vortex dynamics microscopically, the vortex 
velocity- velocity correlation is not calculable directly, be- 
cause the effective vortex Hamiltonian is unknown and is 
precisely what we are looking for. We are forced to aban- 
don the usual approaches of the Nakano-Kubo type, and 
to tackle the problem from the difficult side. Neverthe- 
less, in the DC limit, the transverse force on a moving 
vortex can be calculated from the force-force correlation 
function, in analogy to a DC resistivity formula. This 
limit makes the relaxation time approximation invalid, 
because a significant part of the frequency dependence is 
lost, and the common way of introducing the relaxation 
time approximation by substituting iuj — > icu + 1/t re- 
quires the correct frequency dependence. For example, 
in transport theory, the relaxation time approximation 
is used in an AC conductivity formula, then taking the 
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DC limit subsequently. In addition, the relaxation time 
approximation in a force-force correlation function is al- 
ways erroneous. With an exactly solvable model, we have 
shown that when the relaxation time approximation is 
used in a DC resistivity formula, it leads to results vio- 
lating fluctuation-dissipation theorems. 

Introducing the relaxation time approximation, even 
done correctly, is not a necessary step in obtaining dis- 
sipation. One of the goals of nonequilibrium statistical 
mechanics is to compute various transport coefflcients, 
including the relaxation time, for a given Hamiltonian 
system. In a Hamiltonian system, dissipation appears 
after irrelevant degrees of freedom are integrated out. 
What determine dissipation are quantities like temper- 
ature and strength of impurity potentials, as well as the 
density of state of low-frequency modes of irrelevant de- 
grees of freedom. In the present paper, irrelevant de- 
grees of freedom are the fermionic quasiparticles. When 
those quasiparticle degrees of freedom are eliminated, 
one obtains the vortex friction. The friction formula ob- 
tained hece. follows the one used in dissipative quantum 
dynamics,LJ where it has been explicitly shown that the 
friction obtained by eliminating irrelevant degrees of free- 
dom is equivalent to an evaluation of the random force- 
force correlation function. It also corresponds to the fa- 
miliar Fermi Golden rule for dissipation. We believe that 
a rather detailed study of vortex dynamics based on the 
BCS theory have been presented here, with the key is- 
sue of the sources for the vortex friction. We have shown 
that the vortex friction can come from two contributions: 
At finite temperatures, the finite population of quasipar- 
ticles above and quasiholes below the energy gap give 
rise to a friction which diverges logarithmically at low 
frequency; The nonmagnetic impurities give rise to an 
extra friction which saturates to a value independent of 
the normal-state resistivity in the dirty limit. This core 
state contribution corresponds to the phenomenological 
value obtained in Ref. W. We have also considered the 
effect of coupling to the electromagnetic field and have 
found that it does not change the neutral superfluid con- 
clusions when the superconductor is charge neutral, con- 
sistent with earlier phenomenological treatments. 

Finally, we expect that the method developed here to 
formulate the vortex dynamics in an s-wave superconduc- 
tor will find applications in other systems represented by 
dynamics of collective variables, such as vortex dynamics 
in d-wave superconductors, fission and fusion in atomic 
nuclei, and even the quasiparticle dynamics in quantum 
Hall systems. 
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APPENDIX A: DIVERGENT OVERLAP 
INTEGRALS AND VORTEX FRICTION 

There have been some questions above the implica- 
tions of Eqs. (38) and (39), repeatedly raised by referees 
as well as by others during private discussions, in par- 
ticular on the diverging nature of the overlap integrals 
on the right-hand side of Eqs. (38) and (39), when the 
energy difference between two eigenfunctions vanishes. 
This question may have already been addressed in the 
literature. Nevertheless, we believe it is helpful to give it 
an explicit discussion in the present context. 

We note that Eqs. (38) and (39) are exact conse- 
quences of the fact that the Hamiltonian Ti.Q is the func- 
tion of the parameter xq- To make the connection to 
the scattering problem of quasiparticles scattered off a 
vortex, the thermodynamic limit must be taken first to 
allow the existence of the continuous spectrum. This im- 
plies that it is appropriate to use the Dirac delta function 
normalization for extended states, 



^a = 



Va{x) 



^ifi6 



Vl V2^ 



e' 2 /+,,,,£ (a;) 



(Al) 



the same normalization condition as in Ref. pO. The func- 
tion satisfies the Bogoliubov-de Gennes equation, Eq. 
(30), 



d^ 


1 d 1 (M+i)^ 


-^P 


dr' 


r dr r^ 
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2m 



+ \Air)\f^,^.Eix)=Ef+^^Mx) 



f+.f^.E{x) 
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and 
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Id 1 (/i-i)2 


-P 


dr' 


r dr r'^ 


p 



2m 



+ \A{r)\f+,^,E{x) = Ef^^^,E{x) . 



f-.fj^,Eix) 



(A3) 



Here r — \x~xo\ and k'^+P = kp. Inside the vortex core, 
we may set the energy gap to zero, |A(r)| — 0. There are 
two independent solutions in this region, which we may 
choose to be the following forms: 



24 



A,^,^(x) = i= (^ J ) J^^ , (^^kl + 2m\E\/n'' T 



(A4) 



and 



/2,^,b(x) = -L (^ J J J^_, ^ .^fc2 _ 2m\E\lh' T 



We now consider the left hand side of Eq. (38) with 
an arbitrary small energy difference, 

1 = J d''x¥^{x){W^„no)-^a'ix) ■ (AlO) 

We will show that it can be a finite value (nonzero) . For 
the vanishing small trapping potential Uq, 



(A5) 



(V.o^o) 



Here J„±i are Bessel functions. Away from the zero- 
energy gap region, the corresponding solutions may take 
the forms 







Vx„A 
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Since 



fl4^,E{x) 



V:roA(a;) = -e*''|A(r)|; (x cos 6* + y sin 61) 

-isin^ + ycos0 



-I e 



|A(r) 



(A12) 



J_l ^1+^W^W/E 
and 

f2,t^.E{x) 



J^+i {k+{E)r) , (A6) the integral I may be expressed as 



j_l y i - y/g^ - \A\yE 

V2 [Ji + ^/W^jA^/E 



J^_.{k^{E)r) , (A7) 



where k±{E) = Jkl±2m^E'^ - |A|2/^^- One may 

check that Eqs. (A6) and (A7) give the asymptotically 
exact solutions when r — oo. They are WKB-type solu- 
tions connected to the solutions at r = and r — oo, valid 
under the condition that the energy gap | A| is smooth on 
the scale of 1/kp. Exact solutions may be difficult to find. 
However, for the present purpose of demonstration of the 
diverging overlap integrals they are good enough. The so- 
lutions for a negative energy —E can be constructed by 
using Eq. (37): 



J^x{k+{E)r) (A8) 




and 




J^_i(fc_(i?)r) . (A9) 



The immediate conclusion of the thermodynamic limit 
is that there is an infinite degeneracy for a given en- 
ergy characterized by /x, corresponding to the angular 
momenta of quasiparticles. Those states form the base 
functions for the partial wave analysis of the quasiparti- 
cle scattering, and make the transitions between states 
with the same energy meaningful. 



Jfc^.fcl 
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+ty{±a^yiE,E') + b^^^,{E,E'))] , (A13) 

where x{y) is the unit vector in the x(y) direction, 
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rdr|A|; fl.^Jr)f_^,,^E'{r) 
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and 



b^..f.'{E,E') 



dr\A\ fl^,Eir)f-.^'.E'{r) 



-/-,m,-bW/+,m'.b'W 



(A15) 



The cos6' and sin 9 inside the integral give rise to the 
selection rule for the 

transition elements: J^ '^ d0e~*^^~^ -'* cos ^ = 7r(5p/^^±i 

and J^"" d6'e-*('"-^')^ sin 61 = ±inS^'^f,±i. Only the tran- 
sition between neighboring fi's, that is, /./ — ii ± 1, can 
be nonzero. 

For those nonzero a^^^±i{E,E'), the only possible 
place which may give rise to an infinite value for the in- 
tegral is the region far away from the vortex core. In this 
region, the eigenfunctions are given in the form of Bessel 
functions, Eqs. (A6)-(A(9) (c.f. Eq. (4.10) of Ref. ^, 
which are well behaved. Away from the core |A(r)|J, goes 
to zero rapidly, we conclude that the integral a^^^i {E, E') 
containing |A(r)|^ is finite. It may be instructive to give 
an estimation of afi,^fi,±i{E, E'). For this purpose we con- 
sider /i > and ^' = /i -I- 1 for positive energy states. 
The negative ji and negative energy cases give similar re- 
sults. We will also restrict to the case that E and E' are 
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close to each other, i.e., \E — E'\ < Aqo, and that both 
are sufficiently close to the energy gap, i.e., E ~ Aqo- 
In this case the integral has the largest value, and it in- 
corporates the equal energy limit implied in Eqs. (47) 
and (49). First, we note that since k±{E) k, kp and 
Jf,{k±{E)r) = {k±{E)rY/^i\ for small k{E)r, the Bessel 
function is negligible small ii r < rt — ii/kp and /i is 
large. The integral for region r < r^ is negligible. Sec- 
ond, we also neglect the integral in the region r > ^o, 
because |A(r)|^ is small. The integral now becomes 



a-L, 






(A16) 



Here we have used I to denote the various combinations 
from the solutions, Eqs. (A6)-(A9), specified below. In 
the region r < ^o, |A|^ w Aqo/Co in the above integral. 
For k+{E),k+{E') {k_{E),kJ{E') can be considered in 
the same manner.), following Eqs. (A16) and (A6) we 
have 
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Using the asymptotic form of Bessel function, Ju{z) = 
\jlpKz cos(z— j^7r/2— 7r/4), and approximating the factor 
in the square bracket by 2, we find 
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Here Afc = k+{E) — k+{E'). In reaching Eq. 
(A18) we have dropped a smaller contribution from 



J^^ dr cos[{k+{E) + k+{E'))r - (2^ + 3)tt/2], because 
k+{E)rt = /i > 1, and have also used the fact 
that \k+{E) - fc+(S')l?o < 1- Equation (18) gives 



-'a',m+ 



^{E,E')\<A^/kp. 



For k+{E),k-{E'), following Eqs. (A16), (A6) and 
(A7) we have 
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Please note the difference between Eqs. (A17) and (A19) 
in the indices of the Bessel functions. Since k+{E) — 
k-{E') « 1/Coi and again approximating the factor in 
the square bracket by 2, we find that 
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Co 



A_ 



drcos[{k+{E) - k^{E'))r] 
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Now we consider the phase integral part of I, the inte- 
gral 6^,^, (S, £;'). Forr = |a;-xo| ^ oo, |A| ^ A^o- We 
may ignore the integral in the region r < rt, but not when 
rt > Co- Keeping the leading contribution, h^,^'{E,E') 
may be expressed as 
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The Bessel functions will be replaced by their asymptotic 
forms inside Eq. (A21). In the following we consider 
four cases as done for a^^^'{E,E'). For k+{E),k+{E'), 
following Eqs. (A21) and (A6) we have 
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Since E' —^ E and both are close to the energy gap, 
the factor inside the square bracket is always an order of 
unity, and we approximate by 2. However, we note that 
when E — E' , the term in the square bracket approaches 
to zero when r >> Co- Using the asymptotic form of the 
Bessel functions. 



5l,M 



-Aoo 

A. 



dr- 



1 sin[Afc+-|_ r] 



j.^ 27rfci? r 

sgn{Ak) 7r/2 , \Ak\rt < 1 
, \Ak\rt > 1 



27rfcF I 1/Afcrt 



(A23) 



26 



Here Afc = k+{E) — k+{E'). Again, the contribution 
from J^ drcos[{k+{E) + k+{E'))r - (2/i + 3)7r/2]/r has 
been ignored, because k+{E)rt > 1. Since \k^{E) — 
k+{E')\^o < 1, the condition \k+{E) - k+{E')\rt < 1 
wiU be satisfied if r^ < S.o- 

For k+{E),k^{E'), following Eqs. (A21), (A6), and 
(A7) we have 
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Using a similar procedure for 6i.^ we find 
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We have used k+{E) - k-{E') w l/^o in Eq. (A29). If 
one is concerned about the logarithmic divergence of the 
cosine integral when k+{E) — k^{E') -^ 0, we point out 
that it only occurs when both \E\ and \E'\ are approach- 
ing the energy gap Aqo. In this limit the factor inside 
the square bracket goes to zero linearly, and completely 
removes the logarithmic factor from the cosine integral. 

The conclusion which one may draw from Eqs. (A23) 
and (A25) is that the integral b^,^i±i{E, E') is finite. To- 
gether with what we have obtained for a^.^j,'±i{E, £"), the 
integral I, therefore, the left side of Eq. (38), is finite. 

Using Eq. (38), we have the overlap integral. 
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Since I is finite for the case of /i' = /i± 1 when E^' — Ea — *■ 
0, II diverges as 1/(£'q' — E^)- Because asymptotically 
from the vortex core the wave functions ^q always ap- 
proaches a Bessel function, this diverging behavior may 
be directly deduced from the right-hand side of Eq. (38) 
with the aid of the recurrence relations of the Bessel func- 
tions. The advantage of the demonstration here is that 
the right-hand side of Eq. (49) is finite when two ener- 
gies are exactly equal, without the explicit consideration 
of the diverging behavior of the overlap integral. 

There are two comments worthwhile to make. 

1. The existence of the limit at the left-hand side of 
Eq. (38) in the zero-energy difference indicates that the 
spectral function of Eq. (49) is a smooth function for 
small frequencies, and it may be characterized by a power 
of the frequency. 



2. This limiting behavior also removes the para- 
dox that the frictional coefficient involves inelastic pro- 
cesses, but it may be obtained by calculating the elastic- 
scattering crossection of quasiparticles implied in the 
thermodynamic limiting procedure. 



1 J. Bardeen and M. J. Stephen, Phys. Rev. 140, 1197A 

(1965). 
^ P. Nozieres and W.F. Vinen, Philos. Mag. 14, 667 (1966). 
^ M. Tinkham, Introduction to Superconductivity, 2nd ed., 

McGraw-Hill, New York, 1996. 

* D.J. Thouless, Topological Quantum Numbers in Nonrela- 
tivistic Physics , World Scientific, Singapore, 1998. 

^ P. Ao and D.J. Thouless, Phys. Rev. Lett. 70, 2158 (1993). 
®D.J. Thouless, P. Ao, and Q. Niu, Phys. Rev. Lett. 76, 

3758 (1996). 
■^ G.E. Volovik, Sov. Phys. JETP 77, 435 (1993). 

* E. Simanek, Phys. Lett. A194, 323 (1994). 

®N.B. Kopnin, Physica B210, 267 (1995). This is a clear 
summary of the relaxation time approximation approach 
to vortex dynamics. 

See also 'From microscopic theory to Boltzmann kinetic 
equation: Application to vortex dynamics', ETH preprint, 
1998, by G. Blatter, V.B. Geshkenbein, and N.B. Kopnin. 

^" A. van Otterlo, M.V. Feigel'man, V.B. Geshkenbein, and 
G. Blatter, Phys. Rev. Lett. 75, 3736 (1995). 

" M. Stone, Phys. Rev. B 54, 13222 (1996). 

^^ C. Wexler, Phys. Rev. Lett. 79, 1321 (1997). 

" A.O. Caldeira and A.J. Leggett, Ann. Phys. (NY) 149, 374 
(1983). 

Thouless, Phys. Rev. Lett. 72, 



and D.J. 



■ Q. Niu, P. Ao, 

1706 (1994). 
' D.M. Gaitonde and T.V. 

11951 (1997). 
' J.H. Han, P. Ao, and X. 



Ramakrishnan, Phys. Rev. B56, 



mat/9903125) 



M. Zhu, to be published (cond 



H.B. Callen, Thermodynamics, 2nd ed., Wiley, New York, 

1985. 

M.R. Geller, C. Wexler, and D.J. Thouless, Phys. Rev. 

B57, 8119 (1998). 

R. Kubo, M. Toda, and N. Hashitsume, Statistical Physics 

II, 2nd ed.. Springer- Verlag, Berlin, 1985. 

L.E. Ballentine and W.J. Heaney, J. Phys. C: Solid State 

Phys. 7, 1985 (1974). 

M. Huberman and G. V. Chester, Adv. Phys., 24, 489 

(1975). 

R.S. Fishman and G.D. Mahan, Phys. Rev. B39, 2990 

(1989). 

P. Ao, Phys. Rev. Lett. 80, 5025 (1998). 

P. Ao and X.-M. Zhu, Physica C282-287, 367 (1997). 

J.R. Schrieffer, Theory of Superconductivity, revised, 

Addison- Wesley, New York, 1983. 

X.-M. Zhu and P. Ao, Mod. Phys. Lett. B12 , 451 (1998). 

E.G. d'Agliano, P. Kumar, W. Schaich, and H. Suhl, Phys. 

Rev. Bll, 2122 (1975). 



27 



33 



J.M. Tang and D.J. Thouless, Phys. Rev. B58, 14179 

(1998). 

Y.-C. Chen, J. Stat. Phys. 49, 811 (1987). 

J. Bardeen, R. Kiimmel, A. E. Jacobs, and L. Tewordt, 

Phys. Rev. 187, 556 (1969). 

D. Vollhardt and P. Woelfle, The Superfluid Phases of He- 
lium 3, Taylor and Francis, London, 1990. 

L. Bonig and K. Schonhammer, Phys. Rev. B39, 7413 

(1989). 

R.M. Cleary, Phys. Rev. 175, 587 (1968). 
^* S. Hofmann and R. Kiimmel, Phys. Rev. Lett. 70, 1319 

(1993). 
'■^^ T.A. Brody, J. Flores, J.B. French, P.A. Mello, A. Pandey, 

and S.S.M. Wong, Rev. Mod. Phys. 53, 385 (1981). 
^^ P.G. de Gennes, Superconductivity of Metals and Alloys, 

Benjamin, New York, 1966. 
^'^ G.D. Mahan, Many-Particle Physics, 2e, Plenum, New 

York, 1990. 
^* M. Ma, B.L Halperin, and P.A. Lee, Phys. Rev. B34, 3136 

(1986); S. Mabkawa and H. Fukuyama, J. Phys. Soc. Jpn 

51, 1380 (1981). 
^^ J.D. Jackson, Classical Electrodynamics , 2nd ed., Wiley, 

New York, 1975. 
■*" A.L. Fetter, Phys. Rev. 163, 390 (1967); N. Papanicolaou 

and T.N. Tomaras, Phys. Lett. A179, 33 (1993). 
■*^ Y. Aharonov and A. Casher, Phys. Rev. Lett. 53, 319 

(1984). 



M. Stone, e-print: cond-mat/9708017 



P. Ao, J. Supercond. 8, 503 (1995); Chinese J. Phys. 36, 

190 (1998). 

X.-M. Zhu, E. Brandstrom, and B. Sundqvist, Phys. Rev. 

Lett. 78, 122 (1997). 

J.L Gittleman and B. Rosenblum, Phys. Rev. Lett. 16, 734 

(1966). 



28 



